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USES OF EFFECTIVE FIELD THEORY IN LATTICE QCD 

Andreas S. KRONFELD 
Theoretical Physics Department, Fermi National Accelerator Laboratory a 
Batavia, Illinois, USA 

Several physical problems in particle physics, nuclear physics, and astrophysics 
require information from non-perturbativc QCD to gain a full understanding. In 
some cases the most reliable technique for quantitative results is to carry out 
large-scale numerical calculations in lattice gauge theory. As in any numerical 
technique, there are several sources of uncertainty. This chapter explains how 
effective field theories are used to keep them under control and, then, obtain a 
sensible error bar. After a short survey of the numerical technique, we explain 
why effective field theories are necessary and useful. Then four important cases are 
reviewed: Symanzik's effective field theory of lattice spacing effects; heavy-quark 
effective theory as a tool for controlling discretization effects of heavy quarks; chiral 
perturbation theory as a tool for reaching the chiral limit; and a general field theory 
of hadrons for deriving finite volume corrections. 
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1 Introduction 

The idea to use lattice gauge theory to study quantum chromodynamics (QCD) 
was introduced in 1974 in a seminal paper by Kenneth Wilson. 1 It was an 
exciting time for the strong interactions: it had become clear that quarks are 
partons, 2 the Lagrangian for the gauge theory of quarks and gluons had just 
been published, 3 and the discovery of asymptotic freedom was new. 4,5 Soon 
afterwards an experiment at Brookhaven 6 observed a new resonance, the J, in 
proton-Beryllium collisions, and an experiment at SLAC 7 also observed a new 
resonance, the ip, in e + e~ collisions. This bound state of a charmed quark and 
its anti-quark, now called the J/tp, and its excitation spectrum has gone on 
to play a role to similar to that of positronium in quantum electrodynamics. 
Even today, the charmonium system plays an important role in lattice QCD. 

The emergence of QCD as the theory of the strong interactions means that 
high-energy scattering of partons and bound-state properties of hadrons have 
a common explanation. Owing to asymptotic freedom, 4 ' 5 the gauge coupling 
in QCD becomes weaker at short distances, making pcrturbative cross sections 
more accurate at higher energies. The flip side, however, is that the coupling 
becomes stronger at long distances: perturbation theory breaks down and the 
bound-state problem in QCD is intrinsically non-perturbative. Thus, most 
analysis of non-perturbative QCD has a general nature, relying on, for exam- 
ple, symmetries, analyticity, unitarity, and the renormalization group. At its 
most successful, this kind of analysis yields quantitative relationships between 
experimentally measurable quantities. Perhaps the most striking example is to 
obtain parton densities from deeply inelastic scattering, and then to use them 
to predict jet cross sections in pp collisions. 

Despite such successes, there is a need for general purpose tools to calculate 
properties of the hadrons from the QCD Lagrangian, from first principles. One 
would like to see, quantitatively, that QCD can explain both jet cross sections 
and the proton mass, only by adjusting the QCD coupling and the quark 
masses. One would like to gain insight into the mechanisms of confinement, of 
the type provided by detailed calculations. Finally, to understand interactions 
of quarks at the shortest distances, where new phenomena may be at play, it 
is usually necessary to have calculations of certain hadronic matrix elements, 
with controlled, comprehensible uncertainties. 

Lattice gauge theory provides a mathematically well-defined framework 
for non-perturbative QCD. Before Wilson, Wegner 8 had studied a version of 
the Ising model, on a two-dimensional lattice, with a discrete gauge symmetry. 
Wilson realized how to implement the continuous SU(3) gauge symmetry of 
QCD and other gauge theories and, more importantly, that lattice field theory 
provides a non-perturbative definition of the functional integral. To make any 
sense of quantum field theory, an ultraviolet cutoff must be introduced, and the 
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L = AT s a 

Figure 1: Two-dimensional slice of a lattice with spacing a, spatial size L = Nga, and 
temporal extent L4 = ^40. 



key feature of lattice field theory is that it introduces the cutoff at the outset. 
With the lattice cutoff, it is possible to derive rigorously many properties of 
field theory, as one can see, for example, in the textbook of Glimm and Jaffe. 9 

The basic idea is to replace continuous spacetime with a discrete lattice, 
usually hypercubic, as sketched in Fig. 1. The spacing between sites is usually 
denoted a. For simplicity we consider the spatial volume to have length L = 
N$a on each side, and the temporal extent to be L4 = N^a. From a theoretical 
point of view, the lattice and finite volume provide gauge-invariant ultraviolet 
and infrared cutoffs, respectively. Fermion fields ip(x) and ip(x) live on sites x. 
Gauge fields live on links through the variables 

Un(x) = P exp / ds Ap(x + se M ) (f ) 

Jo 

where P denotes path ordering, and e M is a unit vector in the fi direction. 
In quantum field theory, information is obtained from correlation functions, 
which have a functional integral representation. In lattice QCD the correlation 
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functions are 

(Oi • • • O n ) = i f H dU^x) J] #(z)#(z) Oi ■ ■ ■ O n e' s « CD (2) 

where Z is denned so that (1) = 1, and Sqcd = — J2x £qct>{%) is the (lattice) 
QCD action. The Oi are operators for creating the hadrons of interest and 
also terms in the elcctroweak Hamiltonian. With quarks on sites and gluons 
on links, it is possible to devise lattice actions that respect gauge symmetry. 
As in discrete approximations to partial differential equations, derivatives in 
the Lagrangian are replaced with difference operators. A simple Lagrangian, 
introduced by Wilson, 1,10 is 

£w = -Cwg + £wq (3) 
-^Retr \U„[x)U v (x + ae^Ufa + ae v )Ul{x) - l] (4) 

C Wq = -m^{x)^{x) - Y, fa) [ p +^ ia t - P -n D tut] ^ ( 5 ) 

where g% is the bare gauge coupling, and mo the bare quark mass. (The 
literature often uses (3 = 6/g§ and k = (8 + 2m a)~ 1 .) In Eq. (5) 

P ±/1 = f(l± 7 „), (6) 
D+ U ^(x) = a- 1 [U^xMx + aej - ^(x)] , (7) 

D V lat^( X ) = ^ ~ U U X ~ ae vW X ~ ae »)] ■ ( 8 ) 

The lattice clearly breaks spacetime symmetries, and we shall have to confront 
that issue below. But it preserves SU(3) gauge symmetry, and that is the most 
important feature of Wilson's formulation of field theory. 

It is a simple exercise to show that Wilson's Lagrangian reproduces the 
Yang-Mills Lagrangian in the naive continuum limit, a — > with g% and mo 
fixed. The bare couplings are, of course, unphysical, and the real continuum 
limit must be taken with physical quantities (e.g., hadron masses) held fixed. 
Many discretizations have the correct naive continuum limit. As long as the 
lattice Lagrangian is local, 6 they are all expected to share the same quantum 
continuum limit (with physical masses fixed). The argument for such "uni- 
versality" is based on the renormalization group, which indicates that physics 
depends only only the gauge couplings and quark masses. 11 These arguments 
have been tested experimentally in condensed matter systems, and although a 

''Here "local" means that couplings for interactions of fields separated by a distance r fall 
off exponentially with r. 
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rigorous proof has not been achieved, there is a lot of mathematical evidence 
that lattice field theory defines quantum field theory. 9 

The study of QCD, which a theory of the natural world, is less concerned 
with rigorous theorems than it is with practical results. From this point of view, 
the breakthrough of the lattice formulation is that Eq. (2) turns quantum field 
theory into a mathematically well-defined problem in statistical mechanics. 
Condensed matter theorists and mathematical physicists have devised a vari- 
ety of methods for tackling such problems. In addition to weak-coupling per- 
turbation theory, the toolkit includes non-perturbative versions of the renor- 
malization group, strong coupling expansions, and numerical integration of the 
functional integral by Monte Carlo methods. The first two were pursued with 
great vigor in the decade following Wilson's original paper. The strong cou- 
pling limit is especially appealing, because confinement emerges immediately, 
cf. Sec. 8. The last is the most widely used today, especially for problems 
motivated by particle physics, nuclear physics, or astrophysics. Indeed, when 
most physicists speak of lattice QCD, they mean numerical lattice QCD. 

It is difficult to decide what aspects of lattice QCD to cover in a single 
chapter of a handbook of all QCD. There are several textbooks on lattice gauge 
theory, 12-14 with emphasis on lattice QCD. These books, as well as many re- 
view articles 15 " 19 and summer school lecture series, 20-27 cover the foundations 
well. Nevertheless, many colleagues — theorists whose research is in continuum 
QCD and experimenters whose measurements need non-perturbative QCD to 
be interpreted — tell me that papers on lattice calculations are impenetrable. 

One complaint is that the explanations of numerical techniques are written 
in an unfamiliar jargon, which is unfortunate but hard to rectify here. A more 
serious complaint surrounds the uncertainties that arise in moving the ideal- 
ized problem of mathematical physics to a practical problem of computational 
physics. Many non-experts know that these arise from Monte Carlo statistics, 
non-zero lattice spacing, finite spacetime volume, and unphysical values (in the 
computer) of the quark masses. But, nevertheless, the methods to deal with 
them are (evidently) not transparent. This is a shame. Everyone has a feel for 
statistical errors, even without knowing how to write a Monte Carlo program 
to evaluate the right-hand side of Eq. (2). Similarly, the other uncertainties are 
controlled and understood with effective field theories, so a basis for a common 
language should be possible. 

Two examples illustrate why a common understanding of the uncertainties 
is needed. The first comes from flavor physics, where a wide variety of B, D, 
and K decays are studied, to test whether the standard Cabibbo-Kobayashi- 
Maskawa (CKM) mechanism adequately explains flavor and CP violation. The 
quark-level CKM interpretation of many of these processes is obscured by the 
uncertainty in hadronic matrix elements, of the type (f\H\H), where H is a 
strange, charmed, or 6-flavored hadron, and Ji is a term in the electroweak 
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Hamiltonian. H arises from integrating out W, Z, t, and (possibly) other 
more massive particles. Experimental measurements are, or soon will be, very 
precise. Without reliable theoretical calculations, including a transparent error 
analysis, it will be much more difficult tell whether H is solely electroweak in 
origin or, perhaps, has non-standard contributions. Lattice calculations are 
most straightforward when the final state / has leptons and either one hadron 
or none. c Such matrix elements are needed for leptonic, radiative, and semi- 
leptonic decays, as well as for neutral meson mixing. 

A second example, though less often mentioned, is the search for new 
particles at the Tevatron and, later, at the LHC. If a new particle carries 
color, then its decays always contain jets, and observation depends on the 
extra jet production standing out against normal QCD jet production. The 
QCD uncertainty has, in the past, been estimated by comparing various fits 
to the parton densities. A closer look, 28 however, reveals that the parton 
densities are well-constrained over a limited range of x {x: fraction of the proton 
momentum taken by a parton). Uncontrolled uncertainties in predictions of 
cross sections arise because there are only meager constraints on the parton 
densities outside this range. In lattice gauge theory, however, it is possible 
to calculate the moments of the parton densities: they are related via the 
operator-product expansion to matrix elements of local operators. By the 
time the LHC experiments run, it should be possible to obtain the first few 
moments with an uncertainty that is not only small, but also well justified. 29-32 
The direct calculation of moments and the direct measurement over a finite 
interval in x will provide complementary information, making possible signals 
of new physics more persuasive. 

Suppose the comparison of theory and experiment, in either flavor physics 
or highest-energy collisions, leads to a hypothesis of new, non-standard phe- 
nomena. Then the stakes become very high, and the reliability of error bars 
becomes the central concern. It seems, therefore, worthwhile to discuss how to 
control and estimate uncertainties in numerical calculations, and this issue is 
the central theme of this chapter. The origin of statistical errors in numerical 
lattice QCD is reviewed in Sec. 2. Then the origin of systematic uncertainties 
is reviewed in Sec. 3, motivating the main tools for controlling them. In most 
cases the tool is an effective field theory, which allows us to control the extrap- 
olation of artificial, numerical data to the real world, provided the data start 
"near" enough. Most particle physicists are familiar with the logic and utility 
of effective field theories, and know how to judge their range of validity. Thus, 
effective field theories should also provide a common language for experts and 
non-experts to discuss the error bars, without requiring the non-experts to re- 
peat all the steps of the numerical analysis. The effective field theories needed 



c When there are two hadrons in the final state, there are additional complications, cf. 
Sec. 7.2. 
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to analyze the controllable uncertainties are discussed in Sees. 4-7. 

An uncontrolled systematic effect of many lattice calculations has histori- 
cally been the quenched approximation. It is reviewed in Sec. 2. It is hard to 
estimate the associated error, and only in isolated cases can one argue that it 
is a subdominant error, let alone that it is under control. 

Section 8 contains a "top ten" list of the trends and developments in 
lattice QCD, which warrant appreciation. These have been chosen for their 
broad interest, and because they should influence one's thinking about QCD. 

2 Overview of Numerical Techniques 

This section gives a brief overview the numerical methods. These are covered 
in more detail in some of the texts, reviews, and summer schools cited above, 
as well as in a set of lecture notes aimed at experimenters. 33 

The foremost issue is that there are very many variables. Continuum field 
theory has uncountably many degrees of freedom. Field theory on an infinite 
lattice still has an infinite number of degrees of freedom, but the infinity is 
now countable, i.e., it is as infinite as the integers. This makes the products 
over x in Eq. (2) well-defined. For a computer (with finite memory), the 
number of degrees of freedom must be kept finite. To do so, one must also 
introduce a finite spacetime volume. This may seem alarming, but what one 
has done is simply to introduce an ultraviolet cutoff (the lattice) and an infrared 
cutoff (the finite volume). All calculations in QCD, except trivial ones, require 
an ultraviolet cutoff, and many require an infrared cutoff, although physical 
predictions are cutoff independent. In a sense, removal of the cutoffs is the 
subject of Sees. 4 and 7. 

Even with a finite lattice, the number of integration variables is large. 
For QCD on a Aj x A 4 lattice (cf. Fig. 1) there are (4 x 8)7Vj7V 4 variables 
for gluons and (4 x 3)7VjA4 for quarks. If one only demands a volume a 
few times the size of a hadron and also several grid points within a hadron's 
diameter, one already requires at least, say, 10 points along each direction. In 
four-dimensional spacetime this leads to <~ 32 x 10 4 gluonic variables. 

With so many variables, the only feasible methods are based on Monte 
Carlo integration. The basic idea of Monte Carlo integration is simple: gener- 
ate an ensemble of random variables and approximate the integrals in Eq. (2) 
by ensemble averages. Thus, calling all variables cf), 




z=l 




with weights w to be specified below, and Z defined so that (1) = 1. 
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Quarks pose special problems, principally because, to implement Fermi 
statistics, fcrmionic variables are Grassmann numbers. In all cases of interest, 
the quark action can be written 



where a and (3 are multi-indices for (discrete) spacetime, spin and internal 
quantum numbers. The matrix M a fj is some discretization of the Dirac op- 
erator p + m, such as Eq. (5). Note that it depends on the gauge field, but 
one may integrate over the gauge fields after integrating over the quark fields. 
Then, because the quark action is a quadratic form, the integral can be carried 
out exactly: 



Similarly, products V'a'0/3 in the integrand 0\ ■ ■ ■ O n are replaced with quark 
propagators [M -1 ]^ using the familiar rules of Wick contraction. The compu- 
tation of M _1 is demanding, and the computation of det M (or, more precisely, 
changes in det M as the gauge field is changed) is very demanding. 

With the quarks integrated analytically, it is the gluons that are subject to 
the Monte Carlo method. The factor weighting the integrals is now det Me~ Sg , 
where S g is the gluons' action. Both detM and e~ Sg are the exponential of a 
number that scales with the spacetime volume. In Minkowski spacetime the 
exponent is an imaginary number, so there are wild fluctuations for moderate 
changes in the gauge field. On the other hand, in Euclidean spacetime, with an 
imaginary time variable, S g is real. In that case (and assuming det M is non- 
negative) one can devise a Monte Carlo with importance sampling, which means 
that the random number generator creates gauge fields weighted according to 
detMe~ 5g . With importance sampling the weights on the right-hand side of 
Eq. (9) are independent of the fields, so one can set w = 1. Because importance 
sampling is necessary to make lattice QCD numerically tractable, all numerical 
work is done in Euclidean spacetime. 

Importance sampling works well if det M is positive. For pairs of equal- 
mass quarks, this is easy to achieve. With the Wilson action, Eq. (5), 
75M75 = MK Since this is a similarity transformation, M and M' have 
the same physical content. With M for one flavor and for the other (of 
same mass), the fermion determinant is det(M^M), which is obviously non- 
negative. The same argument holds for Neuberger's discretization, 3 which is 
computationally more demanding, but has better chiral symmetry (cf. Sec. 6). 
For the Kogut-Susskind quark action, 35 the matrix Mks is non-negative, but 
each Kogut-Susskind field creates 4 fermion species in the continuum limit. 




(10) 




(11) 
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Thus, most calculations of detM are for 2 or 4 flavors. The physically 
desirable situation with three flavors, with the strange quark's mass different 
from that of two lighter quarks, is difficult to achieve. One way is to cope 
with occasionally negative weights. 36 Some algorithms for generating the gauge 
fields set up a guided random walk with a finite step size e. 37-39 In a widely 
used scheme for introducing the fermions, 39 they can generate weights such as 
(detMtM) 1 / 2 or (detM KS ) 1/4 , which formally give a single flavor. There is 
some evidence that there may be subtleties associated with non-zero e in large 
systems with small quark masses. 40 These potential problems, and also the 
physical interpretation of the fractional powers, could be monitored by looking 
at the pattern of spontaneously broken chiral symmetry, cf. Sec. 6. 

The choice of imaginary time has an important practical advantage. Con- 
sider the two-point correlation function 

C 2 {t) = (0|$ H (t)*k(0)|0>, (12) 

where §h is an operator with the quantum numbers of the hadron of inter- 
est, H. For simplicity, assume t > and take the total spatial momentum to 
vanish. Inserting a complete set of eigenstates of the Hamiltonian between $jy 
and 

c 2 (t) = J2 m H \H n }(H n \& H \oy mH "\ (13) 

n 

where fnn n is the mass of \H n ), the nth radial excitation with iJ's quantum 
numbers. For real t it would be difficult to disentangle all these contributions. 
If, however, t = ix 4 , with X4 real and positive, then one has a sum of damped 
exponentials. For large X4 the lowest-lying state dominates and 

C 2 (x 4 ) = \m„\H)\ 2 e- m ^ + ■ ■ ■ , (14) 

where \H) is the lowest-lying state and run its mass. The omitted terms are 
exponentially suppressed. It is straightforward to test when the first term 
dominates a numerically computed correlation function, and then fit the ex- 
ponential form to obtain the mass. 

This technique for isolating the lowest-lying state is also essential for ob- 
taining hadronic matrix elements. For a transition from one hadron H to 
another H', one must compute the matrix element (H'\Q\H), where Q is the 
operator inducing the transition. In flavor phenomenology, Q is a term in the 
electroweak Hamiltonian; for moments of parton densities, Q is a local opera- 
tor appearing in the operator product expansion of two currents. One uses a 
three-point correlation function 

C H >Q H (X4,y4) = (0|$^(x 4 + 2/ 4 )Q(2/4)^(0)|0), (15) 
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where only the Euclidean times of the operators have been written out. In- 
serting complete sets of states and taking x 4 and y 4 large enough, 

C H >QH(x4,y 4 ) = mw\H'){H'\Q\H){H\^ H \0)e- m »'^- mH y\ (16) 

The amplitudes (0\® H ,\H'} and (H\& B \Q) and the masses uih 1 and to# are 
obtained from two-point correlation functions C2, leaving (H'\Q\H) to be de- 
termined from Ch'qh- To compute amplitudes for a transition from H to the 
vacuum (as in a leptonic decay), one can simply replace &h in C2 with the 
charged current. 

Equations (12)-(16) assumed t > 0, 24 + 1/4 > 1/4 > 0, and L 4 — > 00. 
With finite L4 other time orderings lead to terms with e ~ m ( L 4-t) They are 
straightforward to derive and to incorporate into fits. Thus, these details do 
not alter the basic paradigm for computing masses and matrix elements. 

These methods are conceptually clean and technically feasible for calculat- 
ing masses and hadronic matrix elements with at most one hadron in the final 
state. The procedure for computing correlation functions is as follows. First 
generate an ensemble of lattice gluon fields with the appropriate weight. Next 
form the desired product 0\ ■ ■ ■ O n , with quark variables exactly integrated out 
to form propagators M _1 . Then take the average over the ensemble. Finally, 
fit the Euclidean time dependence of Eqs. (14) and (16). With two hadrons in 
the final state, correlation functions can be obtained in more or less the same 
way, but the interpretation of the energies and amplitudes is more complicated, 
as discussed in Sec. 7.2. 

Within the same ensemble, there are correlations in the statistical fluc- 
tuations of the quantities calculated. Methods, such as bootstrap and jack- 
knife, that propagate correlations through the analysis are well understood and 
widely used. So, these days, statistical errors rarely lead to controversy. As 
discussed in Sec. 3, it is not practical to carry out the Monte Carlo calculations 
at very small lattice spacings or at very small quark masses. To gain control 
over these effects (using effective field theory to guide extrapolations to the 
physical limit) requires small statistical errors on the raw output of the Monte 
Carlo calculation. 

As mentioned above, the computation of the factor detM in Eq. (11) is 
very demanding. The determinant generates sea quarks inside a hadron. It 
is thus tempting to replace det M with 1 and compensate the corresponding 
omission of the sea quarks with shifts in the bare couplings. This approxima- 
tion is most often called the quenched approximation. 41 A more vivid name 
is the valence approximation, 42 which stresses that the valence quarks (and 
gluons) in hadrons are treated fully, and the sea quarks merely modeled. The 
idea is analogous to a dielectric approximation in clcctromagnctism, and it fails 
under similar circumstances. In particular, if one is interested in comparing 
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two quantities that are sensitive to different energy scales, one cannot expect 
the same dielectric shift to suffice. 

It is not easy to estimate quantitatively the effect of quenching. The 
quenched approximation can be cast as the first term in a convergent 
expansion, 43 providing a method to compute the shifts in the couplings, and 
further corrections. The computed shift agrees with the empirical one, but it 
is about as difficult to compute the next term as to restore the fermion deter- 
minant. For some quantities one can estimate the short-distance contribution 
to the quenching shift. Examples include the strong coupling a s the quark 
masses 45-47 , and which is a form factor needed to determine the CKM 

matrix element 14;,. 48 It is the long-distance part which is harder to fathom. 

The quenched approximation is going away. In heavy quark physics the 
CP-PACS 49 ' 50 and MILC 51 collaborations have unquenched calculations of the 
heavy- light decay constants Jb, /s s , ,/zj, and fjj s . Both groups have results 
at several lattice spacings, so they can study the a dependence. Their results 
are about 10-15% higher than the most mature estimates from the quenched 
approximation. In addition, the Rome group has an unquenched calculation of 
the b quark mass, which agrees well with their quenched calculation. 52 There 
are also unquenched calculations of moments of parton densities. 31 

3 Why Effective Field Theories? 

In this section we discuss why it is necessary, as a practical matter, to consider 
effective field theories. The first clue is that the physical problem has many 
scales, so one should think of renormalization-group strategies to tackle them. 
Effective field theory is one of the most powerful such tools. The central idea 
is to introduce a separation scale (in energy units) \i. Effects from distances 
shorter than are lumped into the couplings (or short-distance coefficients) 
of the effective field theory, whereas effects from distances longer than /x -1 arc 
described by operators in the effective field theory. The degrees of freedom in 
the effective field theory are those required to reproduce the singularities of 
thresholds, 53 and so on, of the underlying theory, when only processes of energy 
E < ji arc considered. By demanding that physical results do not depend on 
yu, and by matching the effective to the underlying theory, one can avoid over- 
or under-counting contributions at the interface of "long" and "short." 

In QCD, the energy scale characteristic of non-perturbative gluonic effects 
is called A, and it lies in a range from the asymptotic freedom parameter 
Aj^jg « 250 MeV through the p meson mass m p — 770 MeV to the scale of 
chiral symmetry breaking m 2 K /m s k, 2500 MeV. For concreteness, we shall 
think of A m 750 MeV, allowing leeway of a factor of three where appropriate. 
QCD also has quarks, whose masses range widely. Light quarks are those with 
m q < A or even m q <C A; the mass of the strange quark is about 100 MeV — 
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seven or eight times smaller than A, and the masses of the up and down quarks 
are about 25 times smaller still. Heavy quarks are those with toq > A or even 
rriQ » A; the mass of the bottom quark is 4-5 GeV — roughly six times larger 
than A, and the mass of the top quark is about 40 times larger still. The mass 
of the charmed quark is slightly larger than A, m c sa 1.3 GeV. Whether, or 
under what circumstances, the charmed quark can be treated as heavy is an 
open question. Systems with valence charmed quarks seem to enjoy some of 
the simplifications of analogous systems with valence b quarks. On the other 
hand, a sea of cc pairs could play a role in the high-momentum tail of hadronic 
wave functions, where not much is known. 

In numerical calculations, one also has cutoffs. The lattice has a non-zero 
spacing a, corresponding to a ultraviolet cutoff ~ n/a. The spacetime volume 
has a finite size L 3 L4, corresponding to a infrared cutoff (L 3 L^) X ^. With 
arbitrarily large computer memory and processing power, one could imagine 
taking these scales in the hierarchy 

L _1 < m q < A < m Q < a" 1 , (17) 

with the quark masses adjusted to their physical values. (Because quarks are 
confined, the adjustment is made by tuning one hadron mass for each flavor of 
quark.) In that case one would only need to know a little about cutoff effects: 
just enough to be confident that the cutoffs introduce only small deviations 
from the limits a — > and L —* oo. 

In practice, however, this idealized situation cannot be achieved. Finite 
computer memory and processor power limit both a and L. To get a feel 
for the tradeoffs, one needs only a few simple scaling laws for the computer 
algorithms. First, the amount of memory needed grows as 

memory oc iVf JV 4 = i 3 i 4 /a 4 . (18) 

Increasing either the physical volume L 3 or decreasing the lattice spacing a 
puts great demands on the memory. The large exponents in Eq. (18) are 
unavoidable, because they stem from the fact that we live in 3 -I- 1 dimensions. 
Second, the amount of CPU time needed to create statistically independent 
lattice gauge fields (with L fixed) grows as 

(19) 

The 4 in the exponent arises because the number of variables to process grows 
as a~ 4 . In addition, the update algorithms slow down as a — > 0, because they 
update in a region of size a, but must propagate these changes over physical 
regions of size A -1 to get a statistically independent gauge field. Thus, the 
exponent z > 0, and available algorithms have z around 1 or 2. 54 
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In addition to creating gauge fields, one must compute quark propagators 
in the background gauge field. These propagators are needed for the valence 
quarks in any hadron. The numerical problem is to solve 

MG = S (20) 

for G, given some source S, where M is the discretized Dirac operator. M , 
G, and S are MxN matrices, where N oc NgN^ is number of quark degrees 
of freedom. The matrix M is sparse (meaning that most of the entries van- 
ish), because it is practical only to put the most local interactions into the 
lattice Lagrangian. The CPU time needed to solve for G, even with the best 
algorithms, is 

r q cx (A max /A min ) p ~ min{l/(m 9 <, {L/af} . (21) 

where A max and A m ; n are the largest and smallest eigenvalues of M, and m q 
is the quark mass. The exponent p depends on the algorithm and is typically 
1 or 2. Equation (20) is also needed in algorithms to incorporate sea quarks 
correctly, 37-39 ' 55-57 where the effective exponent p is 2 or 3. At fixed a it is, 
therefore, costly to reduce m q . For suitable boundary conditions, the volume 
term in Eq. (21) can take over when m q < L _1 , but this regime has significant 
finite size effects, so it is not suited to general-purpose hadron phenomenology 
In addition to difficulties with finite-size effects and with the chiral slowdown of 
the algorithms, statistical uncertainties increase as the quark mass decreases. 
The bottom line is that it is not practical to run the computer at masses as 
small as those of the up and down quarks in nature. 

In typical calculations, these days, Ns — 16-32, with N4 the same or a few 
times larger. To balance the infrared and ultraviolet cutoffs, one ends up with 
a -1 <~ 1-4 GeV (so n/a <~ 3-12 GeV) and L ~ 1-4 fm. In typical calculations 
the "light" quarks have mass in the range 0.2-0.4 < m q /m s < 1.2. Thus, in 
practice, the idealized hierarchy (17) becomes 

< m q < A < m b ~ o -1 , (22) 

which is sketched in Fig. 2. Although the various scales are not as well sepa- 
rated as in the idealized hierarchy (17), there is still some separation. Thus, 
one has a chance of using effective field theory to go from sequences of calcu- 
lations roughly in the hierarchy (22) to continuum, infinite- volume QCD with 
the quark mass of the real world. 

The large masses of the bottom and charmed quarks make the numeri- 
cal solution of Eq. (20) relatively easy, but heavy-quark discretization effects 
become a difficult problem, and a topic of some debate. Instead of studying 
mg <~ a -1 directly, some groups set mq < wif, (and, indeed, mq < m c ), 
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Figure 2: Scales in QCD, and in numerical lattice calculations. Physical scales of QCD are 
labeled on the bottom, and are indicated by solid colors. Scales from practical Monte Carlo 
calculations are labeled at the top, and indicated by cross-hatched patches. 



leading to the hierarchy 

L^ 1 < m q < A < toq < a -1 . (23) 

Another approach for heavy quarks is the static approximation, vtlq — > oo. 
Methods for heavy quarks are discussed further in Sec. 5. 

In summary, practical limitations of computers constrain the lattice spac- 
ing, the quark masses, and the box size to the hierarchy (22) or (23) instead of 
the idealized hierarchy (17). Nevertheless, these parameters can all be varied 
over a certain range, providing numerical lattice QCD with one of its most 
important strengths. The paradigm is as follows: the computer generates nu- 
merical data, varying each of a, mq, m q , and L. These data must then be 
analyzed to extract QCD, at least as long as the data start close enough to 
the real world. With sound theoretical guides, this paradigm is practical, and 
allows propagation of errors. 58 In each case, the guide comes from an effective 
field theory. Moreover, one can test the functional form anticipated from an 
effective field theory and then — assuming the test succeeds — the extrapolation 
to the physical limit is justified. Indeed, this paradigm uses limited computer 
power much more effectively than when relying on brute force alone. 

For reference, the most important effective field theories are listed in Ta- 
ble 1. Quark mass dependence relies on methods that are also used in con- 
tinuum QCD: chiral perturbation theory (%PT) for light hadrons, and heavy- 
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Table 1: Effective field theories used in controlling and quantifying systematic uncertainties 
of lattice calculations. (NRQCD is used for bound states of a heavy quark and a heavy 
anti-quark, called quarkonium, where the small parameter is the relative velocity v between 
the constituents.) 

scale EFT small parameter 

a Symanzik effective field theory Aa 
vhq heavy-quark effective theory (HQET) A/2mQ 

non-rclativistic QCD (NRQCD) v 
m q chiral perturbation theory (xPT) (rriK f^irf^) 2 ~ m s /A 

L Liischer effective field theory e~ m7,L , (LA)' 1 



quark effective theory (HQET) and non-relativistic QCD (NRQCD) for heavy 
quark systems. A glance at Fig. 2 shows that in most cases the scales are 
not very far apart. The leading term in each effective theory may not be 
enough to give a good guide to the extrapolation, but it is possible to work 
out higher-order expressions analytically and use them. 

In some cases it makes sense to combine two or more of the tools. As 
mentioned above, on typical lattices the b quark satisfies rriba ~ 1, so an effec- 
tive theory treating both a and m^ 1 as short distance should be useful. This 
idea is developed in Sec. 5. Most formulations of lattice fermions explicitly 
break some of the (chiral) flavor symmetries of QCD. To study the interplay 
of the continuum limit and spontaneous chiral symmetry breaking, one can 
modify the chiral Lagrangian to parametrize the lattice's explicit chiral sym- 
metry breaking. We shall return to this point in Sec. 6. Finally, when masses 
of pseudo-Goldstone bosons (n, K, rj in nature) become small compared to 
the volume, it becomes necessary to study their finite-size effects with chiral 
perturbation theory. This topic is mentioned briefly in Sec. 7. 

4 Lattice Spacing Effects: Symanzik Effective Field Theory 

The most obvious difference between lattice gauge theory and continuum QCD 
is the non-zero lattice spacing. It is often (correctly) said that lattice field 
theories define quantum field theories. But the definition entails taking a 
sequence of lattice theories, with varying lattice spacing, and taking the limit 
a — > with physical quantities held fixed. 11 To interpret computer calculations 
at non-zero a, however, what one really needs is a description of cutoff effects. 
This section discusses such a description, based on an effective field theory 
invented by Symanzik. 59 ' 60 It provides simple semi-quantitative estimates of 
lattice spacing effects. More interestingly, it provides strategies for eliminating 
them, both by parametrically reducing their size, and by giving a framework 
for combining results from several lattice spacings. 
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One can get some feeling for lattice spacing effects by looking at the quark- 
gluon vertex. With Wilson's action for lattice fermions, Eq. (5), the vertex is 

= -9ot a |7 M cos[i(p + p%a] - ism[±(p + p%a]} 

= -got a {^-^a(p + p') ti + 0(a 2 )}. (24) 

In a hadron at rest, the typical momenta are of order A. So, if Aa is around 
0.2, the quark-gluon interaction deviates by about 20% from the continuum. 

Recall, however, that the lattice action is not unique. In 1985, Sheikh- 
oleslami and Wohlert 61 suggested adding another interaction to the Wilson 
action, namely a lattice approximant to iqa ■ Fq, with coefficient cgw/4. Then 
the quark-gluon vertex becomes 

r^(p,p') = -g t a {7 M cos[i(p + p%a] - isin[±(p + ?J% a ] 
+^c sw a^ v cos[ifc^a] sin[fc„a]} 
- -9ot a {^-%a[(p + p')^+c sw ic7^k»}+0(a 2 )} , (25) 

where k = p' — p, a^ u = i["f fJ- ,"f h ,]/2. On the mass shell, an easy application of 
the Gordon identity shows that the two sets of 0(a) terms cancel if csw = 1- 

The analysis of Eqs. (24) and (25) is essentially classical. For quantum field 
theory, one would rather have a formalism that allows for renormalization, 
and also gives a concept of "on shell" that holds for hadrons. A few years 
after Wilson's paper, Symanzik 59 ' 60 introduced a local effective Lagrangian 
for analyzing discretization effects. This work grew out of his earlier studies 
cutoff dependence in field theories, a line of thinking that had led to the Callan- 
Symanzik equation. 62,63 The idea is that lattice gauge theory, at given a, can be 
described by a local effective Lagrangian. Short-distance effects, in particular 
discretization effects, are lumped into short-distance coefficients, whereas long- 
distance physics is described by local effective operators. The idea is analogous 
to the usage of effective field theories to parametrize short-distance phenomena 
whose microscopic dynamics is unknown. 

For the Lagrangian of any lattice field theory, Symanzik says to write 59 

Aat = Aym, (26) 

where the symbol = can be read "has the same on-shell matrix elements as." 
The left-hand side of Eq. (26) is the lattice field theory inside the computer, 
whose output must be analyzed to obtain a continuum result. The lattice 
action can be complicated, with many free parameters. For lattice QCD with 
Wilson fermions 

Aat = Av + £ a dim °- 4 c Out , (27) 
o 
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where £w is the Wilson (gauge and quark) Lagrangian, Eqs. (3)-(5). The 
Wilson action contains a bare gauge coupling <?q and a (dimcnsionlcss) bare 
mass mod. In Eq. (27), the interactions 0\ at are gauge- invariant composite 
operators of lattice gauge and lattice fermion fields. The couplings Co can be 
chosen in the same spirit as cg-w in Eq. (25), and we shall explain how the 
effective field theory provides a systematic framework for making a choice. 

The right-hand side of Eq. (26) is a local effective Lagrangian (LE£) used 
for analyzing the computer output. Its ultraviolet behavior is regulated and 
renormalizcd completely separately from the lattice of the left-hand side. It 
is convenient to think of £sym as a continuum field theory. The LE£ is the 
Lagrangian of the corresponding continuum field theory, plus extra terms to 
describe discretization effects. For lattice QCD 

£s ym = £qcd + £/, (28) 
where £qcd is the renormalized, continuum QCD Lagrangian, 

C QC j i = ^tr[F^F^}-q(p + m)q. (29) 

The renormalized gauge coupling g 2 and renormalized mass m depend on the 
bare gauge coupling g 2 , the bare mass mo, the co in Eq. (27), and the chosen 
renormalization scheme: 

9 2 = 9 2 (9 2 ,m a;c o - 1 ^a), (30) 
m = m Z m (g 2 ,m a;co; (31) 

where /j, is the renormalization point. The continuum limit is taken for fixed g 2 
and m. Lattice artifacts are described by operators of dimension dimO > 4: 

d = J2 adim0 ^ K o(9 2 ,ma;co^a)0 R (fi), (32) 
o 

where a dim °- i K is a short-distance coefficient, written with factors of a so 
that Kq is dimensionless. As with the renormalized couplings, these short- 
distance coefficients depend on all couplings of the lattice action. 

The renormalized operators Or are sensitive to long distances only. In 
particular, they do not depend on the short distance a. Multiplying matrix 
elements of Or with their coefficients, one finds terms of order (pa) dlmC> ~ 4 , 
where p is a typical momentum, and dimO — 4 > 0. For hadrons consisting 
of quarks and gluons, p ~ A. By assumption Aa is small, so one can treat 
the lattice artifacts in Ci as perturbations. To do so, one can pass to the 
interaction picture driven by £qcd, and in this way develops a series, with all 
matrix elements taken in the (continuum) eigenstates of £qcd- 
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In the interaction picture, one can simplify Li by omitting redundant 
interactions. Let us focus on the quark part of the Lagrangian, and consider 
the field redefinitions 

q^q + a dimX e x Xq, q » q + a dimX e x qX, (33) 

where X is an arbitrary gauge- covariant operator, and Sx and Sx are free 
parameters. Equation (33) simply changes the integration variables of the 
functional integral. On-shell matrix elements, which are the integrals them- 
selves, do not change. 64 Since the interaction picture is being driven by Lqcd, 
the mass shell in question is that of QCD, even though we have not yet solved 
for the hadron masses. 

A trivial example is if X is a constant. Then Eq. (33) changes the normal- 
ization of the fields q and q, so one concludes that on-shell matrix elements do 
not depend on how the field is normalized. In other cases, the field redefinition 
acting on £qcd induces higher-dimension terms, like those in £/. Thus, 



r r i \ ^ dim X 

Li i— » Li + a 

X 



e x qX(p + m q )q + e x q(-p + m q )Xq . (34) 



Similarly, the redefinition acting on terms in Li changes terms of even higher 
dimension. In summary, Eq. (33) amounts to changing certain coefficients 
in Ci. Since the changes are arbitrary, the corresponding operators have no 
effect on on-shell matrix elements. When using the effective field theory to 
describe the underlying theory, their coefficients may be set according to con- 
venience, so they are called redundant. They are easy to identify, because they 
vanish by the equations of motion of QCD. 

Let us illustrate with dimension-five operators. There are no pure gauge 
operators, but there are two linearly independent quark operators, namely 

O b = iqa^F^q, (35) 
0' 5 = 2qD 2 q. (36) 

The second of these can be re- written as 



C>5 = 5 + 2qp {p + m)q- 2mqpq. (37) 

These three terms are, in order, the other dimension-five operator, a redundant 
operator, and something proportional to the kinetic term in £qcd- The last 
can be absorbed into the field normalization of q and a redefinition of m. Thus, 
05 suffices to describe all on-shell dimension-five effects: 



d = aK a . F qia^F^q + ■■■, 



(38) 
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where the ellipsis denotes terms of dimension six and higher. 

The vector and axial vector currents can be described along a similar 
lines. Consider, for example, the flavor-changing transition s — > u. The lattice 
currents take the form 



C = juh^s - acyd^ua^^s + J2 a dim0v - 3 c 0v O v , (39) 

Ov 

Kt = fa-flrf, + ac A d^ u i^ s + J2 a d[mOA - 3 c OA 0% (40) 

Oa 

which are general expressions with the right quantum numbers. In the 
Symanzik effective field theory these currents are described by 65 

Kt = Zv 1 ^ ~ aK v d u ua^s + • • • , (41) 
A L = Za 1 ^ 1 + aKAd^m^s + ■■■, (42) 

where the ellipsis denotes operators of dimension four and higher, ZJ 1 and Kj 
are short-distance coefficients, and 

= u*y*s, (43) 
A 11 = M7^7 5 s, (44) 

are the vector and axial vector currents in continuum QCD. Like the short- 
distance coefficients in the effective Lagrangian, Zj and K j are functions of g 2 
and ma, the lattice couplings c<j, and the parameters co, 7 - Further dimension- 
four operators may be omitted from Eqs. (41) and (42), because they are linear 
combinations of those listed and others that vanish by the equations of motion. 

Like the terms of dimension five and higher in £/, the dimension- four 
currents can be treated as perturbations. Thus, the effective field theory says 
that 

(/latl-^fjiiat) - Wli) + aZ A K A d»(f\ui l5 s\i) (45) 
+ aK a . F J d 4 x(f\T0 5 A^\i) +0{a 2 ), 

and similarly for the vector current. The crucial feature of Eq. (45) is that, 
while the states on the left-hand side are eigenstates of lattice gauge theory, 
those on the right-hand side are eigenstates of continuum QCD. 

The Symanzik effective field theory can be justified in to all orders 
in perturbation theory (in the gauge coupling). In the late 1980's Reisz 
proved a power counting theorem 66 that enabled him to set up a version of 
BPHZ (Bogoliubov-Parasiuk-Hcpp-Zimmermann) renormalization tailored to 
the Feynman diagrams of lattice perturbation theory. With mild assumptions 
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on the gluon and quark propagators, he showed that loop integrals with ex- 
ternal momenta {p} and loop momenta {k} can be written 67 



where Ir denotes parts that are absorbed when rcnormalizing the gauge cou- 
pling and the quark masses. The other term Ijj({p}) does not depend on a 
or on the coefficients co of higher-dimension lattice interactions. Thus, after 
renormalization, lattice perturbation theory is universal: 68 the continuum limit 
does not depend on the discretization. The remainder terms can be developed 
further with BPHZ oversubtractions of the loop integrands. In this way one 
can develop any amplitude's renormalized perturbation series, including con- 
tributions suppressed by powers of a, to any order desired. Similarly, one can 
develop the renormalized perturbation series generated by Symanzik's LE£. 
To lend rigor to the Symanzik theory, one could imagine defining the LE£ on 
an infinitcsimally fine lattice, and using Reisz's methods to extract the univer- 
sal part, and also to define the operator insertions of £/. But in practice any 
method of regulating and renormalizing the ultraviolet will do. 

Although it is only fully justified in perturbation theory, the Symanzik 
effective field theory is believed to hold at a non-perturbative level as well. At 
short distances, small instantons make contributions that should be lumped 
into the short-distance coefficients. But they come in with a high power of 
a and are presumably negligible. Long-distance phenomena, such as confine- 
ment, are more mysterious. But this mystery is bundled into the QCD term 
in the LE£, not £7. Indeed, if the Symanzik formalism were to fall apart, 
it would be hard to understand why other short-distance methods of QCD 
work so well to explain measurements of deeply inelastic scattering, jet cross 
sections, or B decays. 

By calculating on-shell matrix elements on both sides of, say, Eq. (45) in 
renormalized perturbation theory, one can read off the short-distance coeffi- 
cients in £j. In perturbation theory, let 



where n > 0: in some cases the coefficients vanish at the tree level. From 
the justification of the LE£ given above, it is natural to use the renormalized 
coupling in this series. An especially transparent choice of g 2 is something 
physical, at least in perturbation theory, for example the coupling that appears 

d Note that Wilson fermions satisfy the assumptions, but staggered fermions do not. 




(46) 



OO 




(47) 
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in scattering of two quarks of different flavor. Then from familiar properties 
of perturbation theory 

l 

Kq (ma; cp; jia) = K ' r ^ (ma; cp) (ln/xa) r . (48) 

The it dependence must arise in the right way to cancel the /x dependence of 
the operators Or in £j. 

In the discussion so far, the Symanzik effective field theory is simply a 
theory of cutoff effects. 59 It applies no matter how the lattice couplings cp are 
chosen. The default is cp = 0: most higher-dimension lattice interactions are 
omitted from the Lagrangian inside the computer. At this level, the Symanzik 
theory yields, with great sophistication, the result that on-shell matrix ele- 
ments suffer lattice artifacts suppressed by powers of a. 

A lattice gauge theory with smaller short-distance coefficients in the LE£ 
would yield better estimates of continuum QCD. This is where the Symanzik 
formalism becomes especially powerful, because one can demand Kp = 0, and 
solve for the lattice couplings cp. w This is called the Symanzik improvement 
program. A key feature is that if a coefficient Kp (or its expansion coeffi- 
cient Kq) is made to vanish for one observable, then the effective field theory 
set up shows that it vanishes for all observables. 

Of course, it is impractical to solve the equations Kq = in any kind 
of generality. At the very least, one must truncate in (scaling) dimension. 
For QCD, asymptotic freedom guarantees that the anomalous dimensions are 
small, so the scaling dimension is not much different from the classical dimen- 
sion. At dimension five there is only one term in £/, cf. Eq. (38), and with the 
Wilson action, K a .p ^ 0. But, following Sheikholeslami and Wohlert, 61 one 
can write down a discretization of O5, 

5 = tixj^G^TP, (49) 

where is a combination of SU(3) U matrices, defined in Eq. (1), such that 
the naive a — > limit gives F^ u . The Sheikholeslami- Wohlert Lagrangian 

£ S w = £ w + iC S wC 5 . (50) 

The Symanzik coefficient K a .p is, of course, very sensitive to the coupling csw- 
At the tree level 

^fW = j(1-csw)- (51) 

Thus, as in Eq. (25), setting csw = 1 in the computer calculations reduces the 
leading lattice artifact. 

In the foregoing discussion, the short-distance coefficients are written as 
functions of ma, to emphasize that the primary goal of the effective field theory 
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is to separate long- and short-distance scales. For light quarks, however, it is 
safe to expand the short-distance coefficients in powers of ma. One can then 
use this expansion and the ordering of C] by dimension to develop asymptotic 
expansions in powers of a. This is the most widely recognized application of 
the Symanzik effective field theory. Indeed, many reviews of the formalism 
skip over the scale-separation aspect and focus entirely on the lattice-spacing 
expansion. 

To eliminate all terms of order a, the improvement program is relatively 
straightforward. One needs only the first two terms in the small a expansion 
of the normalization factors 

Z.,{ma) =Zj[l + bj\{m u + m s )a] + O(mV). (52) 

For heavy quarks, the expansion in rriQa only makes sense if m,Qa <C 1. As 
discussed in Sec. 3, this situation is not easy to attain for the b quark. We shall 
return to further aspects of the Symanzik LE£ for heavy quarks in Sec. 5. But 
for light quarks Eq. (52) is sensible and useful. In the same vein, it is consistent 
to replace K a .p and Kj with their values at m q a = 0. This set of choices is 
called 0(a) improvement. 

At the tree level, for currents given by just the first two terms of Eqs. (39) 
and (40), respectively, one finds Eq. (51) 

Zf = 1, (53) 
6 [ 7 01 = 1, (54) 
Kf = cf. (55) 

The one- loop corrections to the 0{a) terms have been calculated, with the 
improved lattice Lagrangian £sw, for the Lagrangian itself 69 and for the 
currents. 70 ' 71 The results are not especially informative, except 71 

K [ y ] = c v ] + CV0.01225, (56) 
k a ] = c a ] + CfO.00568. (57) 

Thus, one can obtain the desired improvement condition Ky = Ka = by 
setting c v = - ,g 2 C F 0.01225 + 0{g 4 ), c A = - g 2 C F 0.00568 + 0{g 4 ). 

An important development of recent years are methods for computing the 
short-distance coefficients, through 0(a), non-perturbatively. There are two 
key ingredients. The first is chiral symmetry In 1985 it was observed 72 that, as 
ma — ► 0, both Zy and Za can be computed non-perturbatively by imposing the 
chiral Ward identities. For flavor conserving currents, Zy is simply the factor 
that normalizes the flavor charge. Then, because chiral Ward identities relate 
vector and axial vector correlation functions, they fix Za- The second key 
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ingredient is the Symanzik effective field theory, which provides a framework 
for imposing chiral symmetry through O(a). 65,73 

To compute the 0(a) terms, one can exploit the partially conserved axial 
current (PCAC). 73 In continuum QCD the PC AC relation is 



where the pseudoscalar density V = iuj^s. The implication of writing Eq. (58) 
as an operator equation is that it holds for all combinations of initial and final 
states. A short manipulation with the Symanzik effective theory shows that 



for any initial and final states. By choosing three combinations of states, one 
may use one to eliminate (m„ + m s )Zp/ Za, and then use the other two to 
adjust ca (introduced in Eq. (40)) and csw so that the right-hand side vanishes, 
as desired for continuum QCD. Because the renormalized mass drops out, the 
resulting conditions on ca and csw do not depend on the rcnormalization 
scheme of the effective theory In a non-perturbative calculations, however, 
the 0(a 2 ) terms are ever present, so ca and csw are determined only with an 
accuracy of order Act, coming from the 0(a 2 ) matrix elements on the right- 
hand side of Eq. (59). The improvement coefficient of the vector current, cy, 
and the normalization constants Zj and bj are then determined by imposing 
Ward identities. 

The non-perturbative calculation of csw is straightforward, 74 and the value 
obtained is relatively insensitive to details such as the states used in Eq. (59). 
It also agrees well with perturbation theory. 69 On the other hand, the non- 
perturbative calculation of ca is not so straightforward — two different groups 
find marginal agreement, 74,75 and it has been found 76 to depend on the dif- 
ference operator d£ t . These difficulties are, perhaps, to be expected, since 
in Eq. (59) Ka is multiplied with the small quantity d 2 (f\V\i). The non- 
perturbative results also do not agree well with perturbation theory. The sit- 
uation for cy is similarly unsettled. 75 All estimates for ca and cy yield small 
values, but they multiply large matrix elements. For example, 76 in computing 
f„ the small correction is CAam 2 K /m s , and m 2 K /m s = 2.5 GeV. 

In the non-perturbative improvement program, symmetries of continuum 
QCD were imposed to determine improvement coefficients of the lattice action 
and currents. Another example of the interplay of symmetry and Symanzik's 
theory comes in the lattice calculation of the kaon bag parameter B Kl which 
arise in the theory of K°-K° mixing. Bk is defined by 



dpA" - (m s + m u )V = 0, 



(58) 



(f\Z A d^ t A^ t - (m s +m u )Z P P ut \i) = 

aZ A K A d 2 (f\T\i) + 2aK rT . F (f\ua ■ F l5 s\i) + 0(a 2 ), 



(59) 



(K\Q\K) = lm 2 K f 2 K B K 



(60) 
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where the AS = 2 four-quark operator 



g = S7"(l-75)^(l-75K 



(61) 



and the kaon decay constant Jk ~ 160 MeV is defined through 



(0\A»\K) = ip» K 



(62) 



To calculate fx and Bk one starts with Af^ t and a lattice approximant to Q, 



where Z~j and Cjk are short-distance coefficients. The chiral flavor symmetry 
group in lattice gauge theory is usually smaller than SU(n/) x SU(n/), so 
one must allow for mixing between the target operator Q and several other 
operators, indexed by the subscripts i, j, k. Similarly, the leading lattice 

(7) 

spacing effects are described by several operators Q k . As usual, the operators 
on the right-hand side of Eq. (63) are defined in continuum QCD. Thus, in 
analogy with Eq. (45), the effective theory says 



where the states on the left-hand side are eigenstates of the lattice theory, 
whereas the states on the right-hand side are continuum QCD eigenstates. 
For other four-quark operators the effective theory description also contains 
operators of dimension less than six, multiplied by 1/a. Such "power-law 
divergences" are absent here, because all operators must have AS* = 2. 

Equation (64) holds for all formulations of lattice fermions, but it becomes 
interesting for Kogut-Susskind (or staggered) fermions. One Kogut-Susskind 
field produces four flavors in the continuum limit, so the lists of operators 

f s) 

Q\ contain various flavor combinations. The drawback of four flavors is ame- 
liorated by a remnant of exact chiral symmetry as m — > 0, cf. Sec. 6. The 
symmetry requires K a .p to vanish as m s a. On the other hand, the lattice 
artifacts of the operator do not vanish: the short-distance coefficients Cik are 
non-zero. The matrix elements (K\Qj \K) appearing in Eq. (64) do vanish, 
however, by the flavor symmetry of (four-flavor) continuum QCD. 77,23,78 Thus, 
the leading lattice artifacts in the lattice calculation of ^m K f K B 'k is of or- 
der a 2 . A similar argument applies to the matrix elements of the axial vector 





q=d, s 



(64) 
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current giving fx, so the leading lattice artifact in Bk , when calculated with 
staggered quarks, is also of order a 2 . 

In the last few years, there has been spectacular development in the under- 
standing of chiral lattice gauge theories. 79-81 For vector-like gauge theories like 
QCD, these developments have led to methods with either very small violations 
of chiral symmetry 82 or essentially no violations of chiral symmetry. 83 ' 84 ' 34 
When these methods are used, the leading lattice spacing errors are of or- 
der a 2 . Indeed, one of the methods — domain- wall fermions 82 ' 85 ' 86 — has been 
applied to Bk- 87 ' 88 These calculations also seem to have much smaller 0(a 2 ) 
effects than state-of-the-art calculations 89 ' 90 with staggered fermions. 

With any of several methods — 0{a) improved Wilson fermions, staggered 
fermions, domain- wall fermions, overlap fermions 34 — the leading cutoff effect 
is of order a 2 . To reduce them further, one is confronted with many, many 
operators of dimension six in the Lagrangian, of dimension five in the currents, 
etc. It may not be feasible to eliminate all these effects with non-perturbative 
methods — especially in light of the difficulties with ca- On the other hand, it 
is possible to compute these coefficients in perturbation theory. 

The technical details of lattice perturbation theory are more cumbersome 
that in continuum QCD, as seen in the Feynman rule in Eq. (25). But, since 
the ultraviolet cutoff is built in, the problem is well suited to computer algebra. 
One can generate vertices and propagators and combine them into diagrams 
automatically. 91,92 Although the lattice renders the integrals ultraviolet finite, 
difficulties can arise in the infrared. But these effects are universal, essentially 
by Reisz's theorem. Thus, once the infrared has been understood for a simple 
lattice Lagrangian, the improvement interactions merely add algebraic com- 
plexity that a computer program can handle. These automated methods were 
first introduced for pure gauge theory at the one-loop level, 91 but are now 
being extended to QCD at the two-loop level. 92 

A more controversial issue is the viability of perturbation theory in this 
context. Because lattice gauge theory is a theoretically consistent whole, it 
seems at first glance natural to use the bare coupling g 2 as the expansion 
parameter of perturbative series. This choice ignores the fact that one wants to 
connect lattice gauge theory (at moderate lattice spacing) to continuum QCD. 
The set up of the Symanzik effective field theory, especially as underpinned 
by Reisz's work, encourages the use of a renormalized coupling. Furthermore, 
from a practical point of view, the bare coupling in the Wilson gauge action is 
a disastrous choice: if one compares perturbative and Monte Carlo calculations 
of short-distance quantities such as small Wilson loops, the agreement is very 
poor. On the other hand, renormalized perturbation theory usually describes 
short-distance properties well, especially when the renormalized coupling is 
chosen according to the Brodsky-Lepage-Mackenzie (BLM) prescription. 93,94 
For example, the differences between one-loop BLM perturbation theory for 
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the renormalization and improvement of the vector and axial vector currents 
is easily attributed to (as yet uncomputed) two-loop effects. 95 

Even assuming perturbation theory is accurate, one is faced with a formal 
issue. To illustrate it, let us assume that the we have 0(a) improvement 
of Wilson fermions at the one-loop level. Then the leading lattice spacing 
effects (for light hadron masses) are of order a s (l/a) 2 Aa and (Ao) 2 . Formally, 
the former dominates as a — > 0. Figure 3 shows a s (l/a) 2 Aa and (Aa) 2 as a 
function of a for A = 750 MeV and a s (2 GeV) = 0.2. In the range where 
lattice calculations can be done, a > 0.25 GeV -1 , the "sub-dominant" effect 
is an order of magnitude larger than the "dominant" one. Ideally, one would 
have enough data to fit to both contributors. Otherwise, one is faced with 
the choice of introducing a bias, by ignoring a s (l/a) 2 a and fitting to a 2 , or 
introducing an error, by ignoring a 2 and fitting to a s (\/a) 2 a (or a). Figure 3 
suggests that the uncertainty stemming from the second Ansatz is larger that 
the uncertainty stemming from the first. 

In the last few years, Symanzik improvement has been a major focus of 
research in lattice gauge theory. In light of this renaissance of Symanzik's work, 
it is surprising that many papers still report calculations at only one lattice 
spacing. The improvement program is based on a theory of cutoff effects, which 
clearly demonstrates the utility of repeating the calculation at several lattice 
spacings. To see the benefits, let us consider a simple example. Suppose one 
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has 100 (arbitrary) units of CPU available. Instead of spending all 100 units 
on the finest possible lattice spacing do, one could consider coarser lattices of 
spacing a\ — a / \/2 and a 2 = a /\/2- Spending 65, 25, and 10 units of CPU 
at a , oi, and a 2 would, according to Eq. (19), yield comparable statistical 
errors. Compared to putting all 100 units at oo, the statistical error would be 
a bit larger, but only by a factor of 1.25 (= 1/V0.65). But, if all 100 units are 
spent at do, then one has an uncertainty, of order (Aa) 2 say, which requires 
a guess for the appropriate A: 250 MeV or 2.5 GeV? A calculation based on 
three (or more) lattice spacings does not require a guess, because the added 
information is tantamount to a calculation of the discretization effect. The 
slightly larger statistical error seems a small price to pay. 

5 Heavy Quark Effects: Heavy-Quark Effective Theory 

Some of the most interesting applications of numerical lattice gauge theory 
arise in heavy quark physics. The main aim of experimental B physics is to 
study flavor and CP violation precisely enough to test the CKM mechanism. 
To connect the CKM matrix or, indeed, other short-distance mechanisms of 
flavor and CP violation, one is faced with theoretical formulae of the form 

/ measured \ _ ( kinematic \ / short-distance \ / QCD 
y quantity J y factors J y factor J y factor 

where the measured quantity is a (differential) rate, and the kinematic factors 
consist of measurable momenta and (hadron) masses. Here the short-distance 
factor includes wavelengths less than (100 GeV) -1 ; in the Standard Model, it 
consists of well-determined parameters (like the Fermi constant Gf) and the 
less well-determined CKM matrix. The QCD factor, as a rule, boils down to a 
hadronic matrix element. In the case of \V u d\ and \V US \, isospin and SU(3) sym- 
metries bring the QCD under control, but in other cases lattice calculations are 
essential. 96 ' 97 In some B decays heavy-quark symmetry provides some control, 
but hadronic matrix elements still appear in contributions at the 1/toq level, 
which, generically, are 10% effects. 

Measured in lattice units, the bottom and charmed quarks' masses can 
easily be large. Even if Aa ~ 0.1-0.3, so that the Symanzik effective field 
theory works well for gluons and light quarks, then ra^a <~ 1-2 and m c a about 
a third of that. For this reason it is frequently (but incorrectly) stated that 
heavy quarks cannot be directly accommodated by a lattice. This observation 
overlooks the physical fact that the heavy-quark mass scale is far removed from 
the QCD scale and that, consequently, the dynamics of heavy-quark systems 
simplify. Nevertheless, it is fair to say that lattice spacing effects are more 
challenging for heavy quarks than for light quarks. So, in this section, we 
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discuss how to classify, control, and minimize discretization uncertainties of 
heavy quarks. 

The key is to make use of effective field theories for heavy-quark sys- 
tems. Indeed, from the inception of non-relativistic QCD (NRQCD) 98-100 
and heavy-quark effective theory (HQET), 101-103 these effective field theo- 
ries have been used to treat heavy quarks in lattice gauge theory. Indeed, 
the early papers 98 ' 99 ' 101 ' 102 inspired the development of HQET and NRQCD 
with continuum ultraviolet regulators, as methods for understanding heavy- 
light hadrons 104-106 and quarkonia. 107 More recently it has been shown how 
to use the continuum effective field theories to understand the heavy-quark 
discretization effects of Wilson fermions. 108-111 

To make a connection with Sec. 4, let us start with Wilson QCD and 
examine how the Symanzik theory breaks down when ma 1. A key to 
the Symanzik LE£ is that the leading dynamics are those of Cqcd > while the 
corrections Lj are small. When ma 1, this split into large+small no longer 
holds. First, the expansion of short-distance coefficients in small ma is no 
longer admissible. Furthermore, Ci contains terms that scale as ma to some 
power. Consider, in particular, 

oo 4 

Cl = ...+Y / a" -1 Yl K x ] qXj2(-^a) n q + ■■■, (66) 

X n— 3 M— 1 

which describe deviations from Lorentz (or Euclidean) invariance. (They do 
respect hypercubic rotations.) At each n, the term with /x = 4 is not small, 
because (— "f^D^a) 11 <~ (ma)". But one can cull 74D4 from Ci by applying the 
equation of motion, 112 

-74^40 = (7 • D + m)q. (67) 

Repeated application of the equation of motion to yields a n (7 • D + m) n and 
(nested) commutators of D4 and D. The commutators do not lead to the 
heavy-quark mass, but to the gluon field strength and derivatives thereof. 
Thus, all large terms come from expanding a" (7 • D + m) n and collecting 
(ma) n ~ r qX (7 • D) r q into a new coefficient for qX(-y ■ D) r q. If X = 1 and 
r = or 1, the coefficients in £qcd are modified, and the LEzC takes the form 

£ S ym = £gaugc + <? ^74^4 + ' D + + ( 68 ) 

where the coefficients m\ and \/m\jm2 result from coalescing all terms mul- 
tiplying qq and <pf ■ Dq, respectively. The notation is taken from the energy 
of a quark with small momentum p: 

E( p )= mi + jf- + 0{p i ). (69) 
lm 2 
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Below mi is called the rest mass, and mi the kinetic mass. With these rear- 
rangements, the operators in C' T now all yield powers of Aa, not ma, so they 
still can be treated as operator insertions. 

Equation (68) rests on the same foundation as Eq. (28); only the split 
between large and small is different, reflecting the new situation ma -^C 1. For 
terms in Eq. (66) with I / 1, the rearrangement can be absorbed into the 
short-distance coefficients of L\. For the Wilson fermion action, none of these 
coefficients (except mi) is unbounded as ma — > oo. 101 ' 108,109 It is important, 
when developing an improvement program for heavy quarks, not to sacrifice 
this property. In fact, the currents used in the O(a) improvement program 
discussed in Sec. 4 do not work well for heavy quarks. The foregoing analysis 
even suggests a suitable improvement program: the large-mass behavior of 
C'j remains well-behaved if one mimics it and, hence, omits from Eq. (27) 
operators with extra time difference operators. 108 Similarly, the corrections to 
the currents should not have any time derivatives at all. 110,111 

Unfortunately, unless mi = m 2 the LE£ is no longer "QCD plus small 
corrections," because the normalization of the spatial kinetic energy is incorrect 
by the factor \J m\jmi. In free field theory, for Wilson fermions, 

= 1 - \rn\a 2 + ±m 3 a 3 + • • • , (70) 
m 2 6 

with no term of order ma. This feature persists to all orders in perturbation 
theory. 113,110 The deviations from the desired m\jmi = 1 can be sizable. For 
ma = 0.8 the two terms shown are —0.46 and +0.26. Such strong devia- 
tions from continuum QCD remain when the gauge interaction is turned on. 
The small ma expansions of other short-distance quantities, for example the 
normalization factors of the currents, also break down as soon as ma -ft. 1. 

There are three remedies to this problem. First, one can take ma <C I, 
so that the rearrangement discussed above is unnecessary. Second, one can 
introduce another parameter to the lattice fermion action, so that q^^D^q and 
<77 • Dq can be normalized separately. Third, one can note that it is not lattice 
gauge theory that breaks down when ma ^ 1 but the Symanzik effective field 
theory. This leaves open the possibility that other tools can be used to control 
the discretization effects of heavy quarks. 

Let us start by considering ma <C 1. In principle, this is fine, because one 
can expand the short-distance coefficients in ma, and Eq. (68) reverts to "QCD 
plus small corrections" as in Sec. 4. In practice there are serious difficulties. 
As discussed in Sec. 2 it will not be possible to reduce a enough to make 
mt,a <C I for many, many years. Another way to reduce ma is to reduce the 
heavy quark mass. But if m < mi,, one must use the heavy-quark expansion 
to extrapolate back up to m^. 114 > 115 The simultaneous requirements ma <C I 
and A/m«l fight against each other, making it hard, on accessible lattices, 



Lattice QCD 



31 



to control both systematic errors, not to mention the crosstalk between them. 
To obtain ma < 0.5, one must take m < m c , sometimes as small as 500 McV. 
It is not clear whether this regime can be connected back to mj through the 
1/m expansion. To reach larger quark masses, m ~ 1.5 GeV, ma is sometimes 
as large as 0.7-0.8. Then discretization effects of order (ma) 2 and (ma) 3 are 
large. Finally, the (ma) n errors are amplified in the 1/m extrapolation, but no 
one has a solid idea for estimating how much. Concerns of this kind have been 
voiced before; Sommer 116 and Wittig have insisted on taking ma — > before 
carrying out any extrapolation in 1/m. Then lattice-spacing and heavy-quark 
effects are decoupled, and the main drawback of their analyses (of data in the 
literature), is that the heavy quark mass is too small. 

A variation on this theme is to set up a lattice gauge theory with dif- 
ferent temporal and spatial lattice spacings, a t and a s . Such lattices are 
called anisotropic. The hope 118 is that the heavy-quark mass appears in 
short-distance coefficients as mat, but not as ma s . Then one could take 
a t /a s ~ A/m&, expand the short-distance coefficients in mat, and determine 
the improvement coefficients non-perturbatively. 118 Unfortunately, there is no 
proof that ma s does not arise, and, for Klassen's choice of the lattice cou- 
plings, it does. 119 When ma s does appear in coefficients, one cannot take 
a>t/a s ~ A/mft, and another remedy is needed. 6 

The second remedy is to modify the lattice gauge theory so that the tem- 
poral and kinetic terms are separately adjustable. A simple way to this was 
introduced by El-Khadra et al. 108 Then one could adjust the underlying lat- 
tice parameters so that in the LE£, mi = m-2. The adjustment can be made 
non-perturbatively, by forcing the rest mass and the kinetic mass of a hadron 
to be the same. It works well, 121 but has not been widely used in heavy-quark 
phenomenology In this approach the ma dependence of the short-distance 
coefficients is not as simple as for light quarks. But at least it is possible to 
set m = nib, circumventing the heavy-quark extrapolation. 

The third remedy is to set up the calculation so that heavy-quark meth- 
ods (HQET or NRQCD) and lattice gauge theory work together. There 
are two ways to go about this. One is to derive the heavy-quark theory 
a priori in the continuum, and then replace the derivatives with difference 
operators. 99-103 These methods are called lattice NRQCD and lattice HQET/ 
The other is to note 108 ' 109 that Wilson fermions possess the same heavy-quark 
symmetries 124 ' 125 as continuum QCD. Thus, correlation functions computed 

e Even if they do not tame heavy-quark cutoff effects, anisotropic lattices still can be useful 
for reducing statistical uncertainties, by providing more timeslices in the region of Euclidean 
time where one state saturates Eq. (14). 120 

■^In lattice HQET one treats the 1/m corrections as insertions. 103 The statistical errors are 
smaller, however, if one puts the kinetic term into the quark propagator. 122 ' 123 This is still 
called lattice NRQCD, even when HQET counting is used to classify the 1/m expansion for 
heavy-light systems. 
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with lattice gauge theory can be described a posteriori by HQET (or NRQCD), 
with a logic and structure parallel to the heavy-quark theory for continuum 
QCD. 109-111 Heavy quark symmetry persists for all ma, so the HQET descrip- 
tion can be developed for all ma, with minor modifications that are discussed 
below. These ideas give a systematic procedure for matching lattice gauge 
theory to QCD. It is sometimes called the non-relativistic interpretation of 
Wilson fermions, and sometimes called the Fermilab method. 

For CKM phenomenology, hadrons with one heavy quark are of greatest 
interest. The most important scales are A and the heavy quark mass m. In 
HQET, as used to describe continuum QCD, one separates these two scales 
and, then, treats higher dimensional operators as perturbations, to develop 
a systematic expansion in powers of A/m. For quarkonium — bound states of 
a heavy quark and a heavy anti-quark — there are three important scales, m, 
mv, and mv 2 , where v is the relative velocity between the heavy quark and 
heavy anti-quark. When m is large enough to probe the Coulombic part of the 
potential, v <~ a s (m) is small. Each operator in NRQCD must be assigned a 
power of v, and the effective theory is used to develop an expansion in v and 
a s , which are treated as commensurate. 

Here we would like to use HQET and NRQCD to understand lattice gauge 
theory with heavy quarks (and moderate lattice spacings). As long as mQ 3> 
Aqcd, one can write 109 

Aat = £hq, (71) 

which means that the lattice gauge theory inside the computer can be described 
by a heavy quark effective Lagrangian £hq- The philosophy is in some ways 
similar, but in other ways different from, Eq. (26). The similarity is that we 
would like to use a continuum field theory to describe lattice gauge theory, with 
an eye to understanding and controlling discretization effects. The difference is 
that, for a heavy quark, the descriptive field theory is built from heavy-quark 
fields, not from QCD quark fields. The latter describes both quarks and anti- 
quarks. 126 The heavy quark field, on the other hand, satisfies a constraint, so 
it corresponds either to quarks, or anti-quarks, but not both. The arguments 
supporting Eq. (71) are both concrete, studying the large mass limit of lattice 
gauge theory, 108 and abstract, noting (as above) that the degrees of freedom 
and symmetries are right. 109 

HQET and NRQCD share the same effective Lagrangian, 

£ HQ = ^2 C n t ( m Q>9 2 , m Qa;n/m Q )O n (fi), (72) 

n 

where the C n are short-distance coefficients and the operators O n encode the 
long-distance behavior. The operators do not depend on the short distance 
scales 1/toq or a. It is useful to think of them, as with Symanzik's LE£, 
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as being defined with a continuum ultraviolet regulator, and some convenient 
renormalization scheme. Compared to the HQET/NRQCD description of con- 
tinuum QCD, the main difference is that there are two short distances, I/vtiq 
and a. Because the change is at short distance, the short-distance coefficients 
C^ at must be modified: they depend on niQa, the ratio of short-distance scales. 
Let us recall some aspects of heavy-quark theory. One has 

£ HQ = + + £ (2) + • • • . (73) 

For HQET £hqet contains terms of dimension 4 + s; for NRQCD £nrqcd 
contains terms of order v 2s+2 . In the following, we shall use HQET count- 
ing, but the discussion could be repeated in NRQCD, with straightforward 
modifications. The leading, dimension-four term is 

■^hqet = h v (iv ■ D mi)h Vl (74) 
where h v is a heavy-quark field satisfying the constraint 

ii/ih v = h v , h v itf = h v . (75) 



The choice of the velocity v is somewhat arbitrary. If v is close to the heavy 
quark's velocity, then is a good starting point for the heavy-quark expan- 
sion, which treats the higher-dimension operators as small. The most practical 
choice is the containing hadron's velocity. 

The mass term in £(°) is often omitted. By heavy-quark symmetry, it has 
an effect neither on bound-state wave functions nor, consequently, on matrix 
elements. It does affect the mass spectrum, but only additively. Including the 
mass obscures the heavy-quark flavor symmetry, but only slightly. 109 For two 
flavors, let 9 = (mi c — mib)v ■ x; then the generators 

! i ( e l9 \ 2 i ( -ie w \ 3 i f 1 \ 
T =2\e~ ie )' T = 2\ie~ ie J' T = 2 \0 -I ) ' (76) 

satisfying the SU(2) algebra [r d , r e ] = e d ^ e r^ . When the mass term is included, 
higher-dimension operators are constructed with C = — imiv^. 127 To 
describe on-shell matrix elements one may omit operators that vanish by the 
equation of motion, —ivVh v = 0, derived from Eq. (74). Higher-dimension op- 
erators are, therefore, constructed from V\ = D\ and [D^,V U ] = [D»,D U ] = 

The dimension-five interactions are 

4 1 L T =^ at 2 +4 at O B) (77) 
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where C2 and Cg are short-distance coefficients, and 

2 = h v D 2 ± h v , (78) 
O b = h v s a pB al3 h v , (79) 

with s afj = -ia afj /2 and B a P = rftrfiF' iV . In NRQCD, 2 scales as v 2 , and 
must be treating as a leading term: £nr_q CD = ^hqet + ^2^2- In NRQCD, 
Ob scales as v A 1 as do several operators of dimension six and seven, and this 
collection of operators of order u 4 gives the next-to- leading correction. 100 

At dimension six and higher, many operators arise. The dimension-seven 
Lagrangian contains the first term to parametrize the absence of full rotational 
symmetry; 

4qet - • • • + c^\k J2( D ±) 4h - ( 8 °) 

It is helpful to think of this operator as appearing in the description of contin- 
uum QCD too, but with C^i 1 * = enforced by symmetry. 

One can also develop an effective field theory description of the vector and 
axial vector currents. 110,111 The details have been worked out for decays of a 
heavy quark into a light quark, 110 and for decays of a heavy quark into another 
heavy quark, 111 which is useful for b — ► c transitions. 

We are now in a position to discuss uncertainties in practical calculations. 
The target is continuum QCD, which can be described along the lines given 
above, with different short-distance coefficients. The coefficients are 

m c ° at = to, (81) 

CT* - ^ (82) 

^cont ^(A 4 ) /qo-v 
CB - W (83) 

where to is a renormalized quark mass, and z is a non-trivial function of g 2 
with an anomalous dimension. At the tree level, z = 1. In mass independent 
rcnormalization schemes, the renormalized mass that appears in Eqs. (81)-(83) 
is the (perturbative) pole mass. 

The description of lattice gauge theory with HQET is useful for comparing 
and contrasting the lattice-spacing uncertainties arising in the various heavy- 
quark methods. Since the dependence on rriQa is isolated into the coefficients, 
heavy-quark lattice artifacts arise only from the mismatch of the C„ and their 
analogs C£ ont in the description of continuum QCD. For brevity, we shall focus 
on the three most widely used methods, namely the extrapolation method, 
lattice NRQCD, and the Fermilab method. We shall discuss lattice NRQCD 
first, and then turn to the other two, which both use Wilson fermions. 
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For lattice NRQCD and lattice HQET, the Lagrangian is 

£ gaugc -^D+ lat <f + Y,c n O n (84) 



-lat 



where VP is a two-component lattice fermion field, the O n are discretizations of 
the higher-dimension O n , and the c„ are free parameters. The c„ are chosen 
so that 

->lat ' 



er = ^ (85) 

C? = ^, (86) 

and so on to the desired order. Equation (85) identifies what one means by 
renormalized quark mass; it is obtained implicitly, by adjusting a meson's 
kinetic mass to mg (or my). Equation (86) is matched in perturbation theory. 
(The rest mass is ignored, because it does not affect matrix elements or mass 
splittings.) Solving for the lattice couplings c„ one finds, in many cases, power- 
law divergences as a — ► 0." Therefore, lattice NRQCD calculations must keep 
a <~ 1/toq, and discretization errors are reduced by keeping more and more 
terms in C\ at . One must improve the light quark Lagrangian to the same 
order. The restriction on a is not of much practical importance, because the 
computing challenges discussed in Sec. 2 restrict it to the same range anyway. 

The Fermilab method uses the lattice Lagrangian in Eq. (27) and adjusts 
the free parameters of the lattice action according to Eqs. (85) and (86). The 
solution of these conditions gives the lattice couplings co in Eq. (27) as a 
function of rriQa. In practice, these relations are obtained in perturbation 
theory. The key difference to lattice NRQCD is that, as a — ► 0, the conventional 
Symanzik LE£ also applies. Consequently, the short-distance coefficients of the 
Fermilab method satisfy 

lim = (87) 

Moreover, the corrections to the limiting behavior are related to the short- 
distance coefficients in the Symanzik effective field theory. 110 The pattern of 
uncertainties in the Fermilab method depends on rriQa. For rriQa > 1, dis- 
cretization effects follow a pattern similar to NRQCD, whereas for mqa < 1 the 
Symanzik theory also is valid. On general grounds, one expects the crossover 
region to be smooth, and this expectation has been explicitly verified in several 
cases at the one-loop level. 128,113,110,111 

In the extrapolation method, one (artificially) sets rriQa < 1, and assumes 
that Symanzik improvement is adequate. This leaves 

c lat_ c cont_ (mQa) 2 (g8) 
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Table 2: Parametric uncertainties for heavy quarks. For each method (and heavy-light 
system) the nth row gives the relative uncertainty on Q from Q n . No estimates are given for 
the B system in the extrapolation method, because it is not clear how 1/m extrapolations 
amplify the uncertainties. Numerical estimates, in percent, are made taking a s = 0.2, 
A = 500 McV, trie = 1.25 GeV, m b = 4 GeV. For illustration, we have taken a -1 = 2.5 GeV. 
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1.6 


1.6 


1.6 



for non-perturbative 0(a) improvement. Mass splittings suffer from mis- 
matches of order (m,Qa) n A/mQ = (rriQa) n ~ 1 Aa. For matrix elements, one 
must also look at the normalization factor. Typical recent calculations use a 
normalization factor based on Eq. (52), supplemented with an Ansatz to incor- 
porate full tree- level mass dependence from the Fermilab method. 129-132 This 
leaves an uncertainties of order a s (mQa) n . As discussed above, this would be 
fine if rriQa were small enough. For the charmed quark, a preliminary study 
shows that these effects are under control for the spectrum, if one takes the 
continuum limit. 133 

To get a semi-quantitative feel for the uncertainties, let us consider a 
generic quantity Q with heavy quark expansion 

Q = Qo + — + -%. (89) 
m Q m 2 Q 

Table 2 lists the relative uncertainty on Q from each term. The heavy quark 
expansion is useful here, because in all methods the physical heavy-quark ef- 
fects are intertwined with lattice spacing effects. We shall assume that Eq. (89) 
is adequate for charmed hadrons, but the conclusions do not really depend on 
the assumption. For lattice NRQCD and the Fermilab method, we imagine 
that most of the normalization of Qo is non-perturbative, 134 ' 135 ' 48 but the rest 
is available at the one-loop level, and that Qi is normalized at the tree level. 
For the extrapolation method, we imagine full O(a) improvement, and neglect 
the difficulties with ca and cy mentioned in Sec. 4. Table 2 also includes 
numerical estimates, made taking a s = 0.2, A = 500 McV, m c — 1.25 GeV, 
m b = 4 GeV. 

For B physics, lattice NRQCD and the Fermilab method lead to the same 
estimates. With lattice NRQCD, heavy quark propagators require negligible 
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computing; with the Fermilab method, they require more computing, but the 
overhead is negligible compared to generating gauge fields. The main advan- 
tage of the Fermilab method is that it has no restriction on rag a. No esti- 
mate is given for the discretization errors for B physics from the extrapolation 
method. One has to understand how errors of order (ma) n propagate through 
the 1/m expansion. One crude method is to compare different fits; this only 
tests whether the function is smooth in the region where data are available 
and is, thus, an underestimate. 

For D physics, the 0(a) method and the Fermilab method have similar 
uncertainties, because m c a < 1. As soon as m c a < a s , non-perturbative 
improvement becomes valuable. It could be implemented for the Fermilab 
currents without any conceptual difficulty. Lattice NRQCD is rarely used for 
charmed hadrons, because a -1 ~ m c and then perturbation theory in a s , 
needed to improve the NRQCD Lagrangian to interesting accuracy, does not 
converge well. 

To reduce the uncertainty from heavy-quark discretization effects over the 
short term, one probably requires (automated 92 ) perturbation theory In lat- 
tice NRQCD and the Fermilab method, one more loop would reduce each 
uncertainty by about one fifth. Attacking B physics in this way, one could 
direct increases in computing at removing the quenched approximation. In 
the extrapolation method, the program to reduce uncertainties is to reduce a 
until there is a window with rriQa <C 1 and A/rriQ <C 1 simultaneously. Un- 
fortunately, this choice postpones unquenched calculations for at least another 
generation. 

6 Light Quark Effects: Chiral Perturbation Theory 

Light quarks bring physical scales m q <C A into QCD. The numerical algo- 
rithms for computing the quark propagator slow down for light quark masses, 
as explained in Sec. 3. This problem makes it impractical, as a rule, to carry 
out numerical calculations with masses as small as those of the up and down 
quarks. To reach the physical region, the Monte Carlo is run at a sequence 
of light quark masses, say in the range 0.2m s < m q < m s , and masses and 
matrix elements are extrapolated to down to m,d and m u . 

Because of the artificially large light quark masses, hadrons in the com- 
puter carry a cloud of light pseudoscalar mesons, all with mass squared 
nipj < 2m 2 K . In nature, the cloud contains pions, kaons, and r\ mesons, 
with a wide range of masses, m 2 -C m 2 K w \rv? rj . The relation between the two 
situations can be understood quantitatively, by studying the interaction of the 
pseudoscalars with other hadrons. If the numerical data are close enough to 
the chiral limit, then the machinery of chiral perturbation theory 136 ' 137 (%PT) 
can be applied. Numerical data can be tested against the leading (or next-to- 
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leading or ncxt-to-next-to-leading) order prediction of %PT. If the data verify 
%PT, one has a sound guide to extrapolate in the quark mass. 

Chiral perturbation theory is developed by introducing an effective La- 
grangian built of pion fields and, where appropriate, other hadrons. In this 
case, the long distances are the inverse light quark masses, whereas the short 
distance is the QCD scale A. The separation scale \x is, therefore, usually taken 
to be around 1 GeV. Interactions in the chiral Lagrangian can be classified by 
the number of derivatives and powers of the quark mass. The momenta of 
pseudoscalars and the virtuality of other hadrons is also assumed to be small, 
q 2 <~ m 2 K C A 2 . A concrete scale, less fuzzy than A, that arises in the expan- 
sion is A-kJk = 1.65 GeV. 

This chapter is not the place for a full review of xPT. Instead, we shall 
take a specific example — the calculation of the B meson decay constant — 
to illustrate how xPT for lattice gauge theory differs from phenomenological 
applications to low-energy hadron physics. In particular, we will contrast usual 
%PT with two situations that arise in lattice calculations. The first is quenched 
xPT, which is used to describe problems in the quenched approximation. The 
second is partially quenched xPT, which applies to computer calculations in 
which the valence quarks have a different mass than the sea quarks in loops. 
Finally we discuss how to modify the chiral Lagrangian to take into account the 
fact that lattice QCD usually has less chiral symmetry than continuum QCD. 

The classic example of a chiral extrapolation is shown in Fig. 4. The 
pseudoscalar mass is related to the quark mass, in leading order, by 

m 2 PS = (m q + rriq)B, (90) 

where B, from the point of view of chiral perturbation theory, is an unknown 
"low energy constant." From the point of view of lattice QCD, it is calculable. 
One simply computes nips for a variety of (light) quark and anti-quark masses 
and, if the behavior in Eq. (90) is verified, fits to obtain B. Figure 4 shows that 
Eq. (90) is a good description out to remarkably large quark masses. (This 
may be a fluke of the quenched approximation.) Indeed, the lattice calculation 
of m = \{vn u + m d ) is nothing but m = m 2 /B, with = 140 MeV from 
experiment and B from Fig. 4. 

The mass of the strange quark is determined in a similar way. One can 
set Tfiq = rriq and identify m s as the quark mass giving m 2 PS = 2m 2 K . Al- 
ternatively on can extrapolate in the anti-quark mass, until mq = m <C m s 
and identify m s when m 2 PS = {m 2 K+ + m 2 KO )/2. As a rule, in lattice QCD we 
neglect electromagnetic effects and isospin violation — other uncertainties, even 
statistical errors post-extrapolation — are larger. Note that methods have been 
developed to include these effects, when necessary. 139 

Of course, the renormalizcd quark mass (and, hence, B) depends on 
the renormalization scheme. It is important to take this scheme depen- 
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Figure 4: Plot of the squared pseudoscalar mass m 2 ps vs. quark mass m q , with numerical 
data in the quenched approximation. 138 



dence into account before quoting a numerical value. The two most widely 
used conventions are the renormalization-group invariant mass 1 and the MS 
mass fh q {2 GeV). In mass ratios the scheme dependence cancels, and for this 
reason it is convenient to refer to light quark masses in units of the strange 
quark mass, as defined in the previous paragraph. 

For other masses and matrix elements it makes more sense to eliminate the 
scheme-dependent quark masses in favor of pseudoscalar meson masses. The 
meson masses are what one computes numerically, and they are what appear 
in xPT. For example, let is consider fs q , where q is the flavor of the light 
quark in the B meson. We shall neglect 1 /mi, corrections and write 



[1 + A/b,] 



(91) 



where $ is independent of both heavy and light quark masses, and A/s 
denotes the (one-loop) contribution of the light meson cloud. 

In QCD, the one-loop correction to the decay constants are (neglecting 
isospin breaking) 141,142 



A/ Bs = - 



l + 3, 9 2 

(47T/) 2 



[m 2 K ln(m 2 K /fi 2 ) + ln(m 2 /^ 2 )] 



+ d(fi)(m 2 K + \m\) + c 2 (/i)(r 



K 



X). 



(92) 
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AfB « = ~JX^h ln ( m '/^ 2 ) + W K \n{m 2 K /^ 2 ) + ^m 2 ln(m 2 /^ 2 )] 



where / and g are (the chiral limit of) the pion decay constant and B-B*-tt 
coupling The "low-energy" constants c» encode QCD dynamics from distances 
shorter than yu -1 , whereas the logarithms are long-distance properties con- 
strained by chiral symmetry. The dependence on [i cancels in the total. 

These formulae illustrate how to derive useful formulae from xPT. The 
light quark masses in numerical lattice calculations are all about the same size, 
say m q > 0.2m s , so it is not helpful to neglect m 2 relative to m 2 K . Because each 
quark mass is an adjustable parameter in lattice calculations, general formulae 
with non-degenerate quark masses arc needed. If such formulae are available, 
the combination of lattice calculations and xPT is powerful. Chiral symmetry 
constrains the coefficient of the logarithmic terms, tying them to quantities 
that can be calculated by other means. This is a key, because it is difficult, 
from fitting, to distinguish m 2 In m 2 from a polynomial, unless the range of m 2 
is very wide. With the "chiral logs" constrained, however, a straightforward 
fit the determines the low-energy constants. 

Quenching changes Eqs. (92) and (93) drastically. Figure 5 shows the 
quark flow of several virtual processes that take place in a meson. The 
quenched approximation includes Fig. 5(b), but not (a) or (c). As a con- 
sequence, some pion loops are omitted, and rj loops are mistreated. 144 In a 
partially quenched calculation, Fig. 5(a) and (c) are restored, but the masses of 
quarks in loops are not the same as those on the valence lines. In the following, 
we shall examine how the omission or modification of these processes changes 
the quark mass dependence of Jb- 

Recall that the quenched (or valence) approximation replaces the fermion 
determinant in Eq. (11) with 1. This replacement omits closed fermion loops 
(and absorbs some of the omission into the bare couplings) . Instead of omitting 
the loops, one could imagine canceling them with bosonic loops with the same 
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Figure 5: Quark line configurations that lead to meson loops. 
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action. Written as functional integral 145 
Det(p + m) 



Det(p + to) 



= J VqVqDqVqe-^ +m ^ q -^ +rh ^ (94) 



where q is a bosonic field, called scalar quarks or ghost quarks. We want to 
understand partially quenched calculations with n / flavors of sea quarks and 
n v flavors of valence quarks. The mass matrices are 

to = diag(mi,mi, . . .,m nf ,M 1 ,M 2 , . . . ,M n J, (95) 
to = diag(Mi,M 2 ,...,M n J, (96) 

so the ghosts cancel the determinants of the valence quarks. When studying 
mesons, one wants n v — 2; baryons, n v — 3. The quenched approximation is 
then the special case nj = 0. 

The advantage of Eq. (94) is that it immediately suggests the hadronic 
content of the quenched approximation. 145 In addition to the usual qq mesons, 
one has bound states of ghosts qq and bound states of quarks and ghosts qq 
and qq. States with an odd number of ghosts have negative metric and some 
spooky consequences. 

If to and m are both small compared to A, there is an approximate chiral 
symmetry. It is not so straightforward to identify the group, because of dif- 
ferences between bosonic and fermionic integrals, and one must insist that the 
determinants cancel. 146 ' 147 It turns out that, for rif > 1, that the symmetry 
group has the same number of generators as 

Gpq = SU(n/ + n v \n v ) x SU(n/ + n v \n v ) x U(l), (97) 

and the Ward identities derived with this "intuitively obvious" group are the 
same as for the correct group. Here SU(ri2|ni) denotes a graded Lie group, a 
mathematical beast that also appears in supersymmetry. The SU groups are 
spontaneously broken to their diagonal subgroup, just as in QCD. 

For the quenched approximation, rif = 0, the result is slightly different 147 

g q = [sukk) x sukk)] x u(i), (98) 

where 1X1 denotes a semi-direct product. The semi-direct product M may be 
less familiar that the direct product x . Let A and H be subgroups of G, and 
suppose A is normal. (A is normal if, for all a e A, gag^ 1 E A for all g e G.) 
If, furthermore, A n H = {1}, then the semi-direct product A XI H consists 
of all ah, a e A, h e H. This seemingly technical issue has a crucial physical 
implication: in quenched xPT the flavor singlet mesons (r/ in QCD), as well 
as singlet ghost-quark and ghost-antighost mesons, do not decouple. 
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Consequently, for quenched xPT, one must construct the effective La- 
grangian for pseudoscalar mesons with a formalism that includes the r/'. This 
can be achieved 148 along lines similar to normal xPT. 137 The main difference is 
that the r\' is not a normal particle. Its propagator has a double-pole structure, 
because Fig. 5(c) is absent. On the other hand, although it was not obvious 
until the technical details were sorted out, the flavor singlets do decouple once 
there are closed quark loops to restore Fig. 5(c). 147 

Let us now quote the results from quenched and partially quenched x?T 
for A/b- In the quenched approximation 149 ' 142 



Af- Q = — 

J bv ( 4 ^J)2 



3s 2 . 2 
m + jm vl 



\n{ml v /n 



+ c?(/x)m^ + c$(v)m 2 vv , (99) 

where 7 is a combination of low energy constants of the quenched theory, and 
the mass of the light pseudoscalar made from valence quark v. The 
constant m 2 , is the residue of the rj'-like double pole. Finally, and are 
counter-terms, whose dependence on \i renders A/^ independent of fi. In 
general, the low-energy constants of the quenched approximation cannot be 
related in any rigorous way to those of QCD. 

The terms proportional to m 2 diverge in the chiral limit m\ v — > 0. This 
behavior shows in detail that the quenched approximation cannot be forced to 
account for a wide range of scales. The common practice to obtain fg in the 
quenched approximation does not follow from a fit to Eq. (99). Instead one 
assumes that there is little error from quenching when m v ~ m s , and then 
extrapolates linearly in rr? vv . Thus, common practice circumvents the clearly 
unphysical quenched chiral log, but it must be admitted that the uncertainties 
from quenching and chiral logs become intertwined. 

In the partially quenched case, the chiral logs behave better. If all n/ sea 
quarks have the same mass 142 



AfPQ = l + 3g 2 

Jbv - ( 47r/ )2 



-—J-ln(m vf /n ) + JJ 2 ^ Hm vv /fi ) 



+ ci{n)m ff + c 2 {n)m vv , (100) 

which is qualitatively like Eqs. (92) and (93). A simple argument relates the 
partially quenched low-energy constants to QCD. 150 Consider any matrix ele- 
ment to be a function of n/ + n v sea and valence masses. If the valence quark 
masses are set equal to the sea quark masses, then one recovers QCD. Thus, 
QCD lies on a hyperplane in the space of all masses. If all masses are sent to 
zero together, there is a unique chiral limit. Thus, the low energy constants of 
the partially quenched theory and QCD are the same. 
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The foregoing analysis discussed chiral symmetry as it appears in the con- 
tinuum limit of QCD, and its quenched and partially quenched theories. The 
two most widely used formulations of lattice fermions — Wilson fermions and 
staggered fermions — have a smaller symmetry group. Inspection of Eq. (5) 
shows that two sources break SU(n^) x SU(n/), namely the mass term and the 
part of the difference operator proportional to Dirac matrix 1. Wilson added 
the extra terms (or, equivalently, chose such a peculiar difference operator) to 
circumvent the so-called doubling problem. 

The first action that Wilson tried (and, independently, Jan Smit) was 

£ Nq = -m 'tp(x)ip(x) (101) 
- — ^ ${x)l» [U^(x)tP(x + ae^) - Ufa - ae^(x - ae^)] 

If m = this action is chirally symmetric. Unfortunately, the propagator has 
many poles 

S(p) = — -, (102) 

1 L p 7 M sm(p^a) 

where components of the momentum four- vector lie in the range — tt/cl <p^< 
n/a. There is a pole at all 2 d combinations of p^a — 0, ir. These correspond 
to 2 d species of physical particles: vacuum polarization, for example, gets a 
factor 2 d . Moreover, the axial symmetries are exact: the different species 
have a pattern of axial charges such that the anomaly cancels. To avoid these 
problems, Wilson decided it was less drastic to break the axial symmetries 
explicitly, leaving SU(n/) as the flavor group. 

The staggered formulation of lattice fermions starts with Eq. (101), but 
notes that similarity transformation 151 ' 152 

ip{x) i — ► T(x)ip(x) (103) 
ip(x) i — ► i>{x)T\x) (104) 
T{x) - (105) 

diagonalizes all Dirac matrices. One then has a sum over four equivalent terms. 
One can three of them, leaving 2 4 /4 = 4 species in the continuum limit. The 
remaining species have a remnant axial U(l) symmetry 152 On the other hand, 
the flavor group is smaller than SU(4), and the notion of flavor is tied up with 
spacetime symmetries — hence the name "staggered." The SO (4) Euclidean 
invariance is broken by the lattice down to the semi-direct product T 4 1X1 SW4, 
where I^ is the Clifford group of 4 x 4 Dirac matrices, and SW4 is symmetry 
group of a hypercube. 153 ' 154 The factor Y± is interpreted as a flavor group, and 
SW4 as a spacetime symmetry. 
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Table 3: Pattern of chiral symmetry breaking for various formulations of lattice fcrmions. 



formulation 


G 


-» i? 


Wilson 


SU(n/) 


- SU(n/) 


staggered 


u(i) x r 4 


^r 4 


Ginsparg- Wilson 


SU(n/) x SU(n/) 


-» SU(n/) 


continuum QCD 


SU( n/ ) x SU( n/ ) 


-» SU(n/) 



Before coming back to the implications of the reduced chiral symmetry, let 
us review the nature of the problem. The Nielsen-Ninomiya theorem states that 
it is impossible to maintain chirally symmetry and avoid species doubling with 
ultralocal interactions. 155 ' 156 Here uZtralocal means that the lattice couplings 
vanish if the fields are separated by more than a few lattice spacings. The way 
out, only recently appreciated, is to forgo ultralocality in favor of locality, 
which only requires that the couplings fall off exponentially with separation. 
Then, if the lattice Dirac operator p satisfies 157 

j 5 p + Pl5 = ap~/ 5 p, (106) 

which is called the Ginsparg- Wilson relation, then correlation functions are 
chirally symmetric. There are two known solutions to Eq. (106), one based on 
renormalization group ideas related to the original derivation, 84 and the other 
related to the Narayanan-Neuberger formulation of chiral lattice fermions. 34 
Both are local, so they are the basis of a well-behaved field theory, but not 
ultralocal, so they avoid the hypothesis of the Nielsen-Ninomiya theorem. The 
solution of Eq. (20) for these formulations is much more computationally de- 
manding than for Wilson or staggered fermions. Consequently, they are only 
beginning to make their way into numerical calculations. 

Table 3 summarizes the flavor symmetry group for lattice fermions. For 
staggered fermions, the flavor group comes in a semi-direct product with the 
symmetry group of the hypercube, SW4; for the others it comes in a direct 
product. The short-distance coefficients of the chiral Lagrangian must be mod- 
ified to depend on aA, the ratio of two short distances. Because the Symanzik 
effective field theory shows that the violations of chiral symmetry come in 
through higher dimension operators, an expansion in aA presents no difficulty, 
except to introduce new low-energy constants. The details depend on how the 
symmetry is broken. They have been worked out for Wilson quarks 158 ' 159 and 
staggered quarks. 160 

Now that many groups are acquiring the computer resources needed for 
unquenched QCD, the topics discussed in this section will gain in importance. 
Checks of chiral behavior will be important not only for extracting physics from 
the numerical calculations, but also for testing the algorithms. For example, 
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Sec. 2 noted that odd numbers of flavors are sometimes incorporated into the 
Monte Carlo by generating a weight proportional to a fractional power of a 
determinant. For example, with staggered fermions one has (det Mks)™^ 4 , 
because M K s corresponds to four flavors. It is conjectured, 161 that the appro- 
priate graded group for the continuum limit of this theory is SU(4|4 — rif), like 
a partially quenched theory with n v = 4 — n/. This conjecture, if true, has 
important consequences for testing the physical content of these algorithms, 
which are otherwise a bit mysterious. 

7 Finite Spacetime Volume Effects 

In this section we address systematic effects stemming from the infrared cutoff 
imposed by the finite spacetime volume. In Euclidean space, the temporal 
direction starts out on the same footing as the spatial directions. Even so, 
the physical interpretation of finite temporal extent, on the one hand, and 
finite spatial volume, on the other, is different. It turns out that both kinds of 
effects are valuable, providing tools to extend the range of problems to which 
Euclidean, numerical lattice QCD can be applied. 

In Sec. 7.1, we shall address the implications of the finite temporal extent, 
and the relation to QCD thermodynamics. In Sec. 7.2, we turn to the effects 
of the finite spatial volume. Several issues arise here. Generic corrections to 
masses and 1 — » and 1 — » 1 matrix elements are suppressed by a factor e~^ L , 
where fi is a mass related to the mass of the lightest hadron (assumed massive). 
The effects are larger, suppressed by powers of L instead of an exponential, 
when pseudo-Goldstone bosons satisfy Lm, < 1 . The most intriguing effect of 
finite volume is on scattering states, whose allowed energies are connected in 
a model-independent way to final-state phase shifts. 

7.1 Finite temporal extent (non-zero temperature) 

Recall that numerical lattice calculations are formulated in Euclidean space, 
namely with imaginary time. In the computer, the extent of the imaginary 
time is finite, L4 = N 4 a. In Sec. 2, starting with Eq. (12), we blithely assumed 
a correspondence between correlation functions and vacuum expectation values 
of time-ordered products, viz., 

\ f n^wn^ww 0i(ti)0 2 (t2)---0n(tn) e ~ s ^ 

X,fJ, X 

(0|T O l (t l )0 2 {t 2 )---O n (t n ) |0). (107) 

When L4 — N 4 a is finite, however, this correspondence is not exact. 

The essential features are easily seen for continuous time. Then the Eu- 
clidean action is simply related to a Lagrangian (of a system with several 
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degrees of freedom at each lattice site) . From this Lagrangian one can straight- 
forwardly derive a Hamiltonian and set up the equivalence between the path- 
integral and canonical formulations of quantum mechanics. To be explicit, 
consider the case with only one operator 0{t) in Eq. (107). Also let us start 
with fixed boundary conditions, which means that fields at x 4 = and X4 = L 4 
are not integrated over. Then, 

(O(t)) = |(/|e-^-*)6e-^K), (108) 
Z = (f\e~ AL *\i), (109) 

where \i) and (/| are the initial and final states implied by the boundary condi- 
tions, H is the Hamiltonian, and O is the Hilbert-space operator corresponding 
to the function (of quark and gluon fields) O on the left-hand side. Inserting 
complete sets of eigcnstatcs of H into Eq. (108) one sees 

£n e ~ E " L *{f\n){n\i) 

where the sums range over all eigenstates of the Hamiltonian, and E n is the 
energy of state \n). For large L 4 the ground state — or vacuum — dominates the 
sum in the denominator: 

Z->e-^</|0)<0|i). (Ill) 
Similarly, if t and £4 — t are large 

(O(t)} = (0|O|0) + O ( e - {El - Eo) mi «( L 4-M)) ; (112 ) 

with the vacuum energy and overlap factors canceling out of the ratio in 
Eq. (110). One can repeat this analysis for operators at several times and, 
in this way, derive Eq. (107). 

Although fixed boundary conditions are sometimes useful, 162 ' 163 the con- 
tamination from excited states is unnecessarily large. This is especially so in 
two- and three-point functions, where one wants several time separations to be 
large, cf. Eqs. (12)— (16) . To reduce the contamination term in Eq. (112) it is 
helpful to choose boundary conditions so that |/) = e~"^|i), and then to sum 
over i. Then Eq. (110) becomes 

«*» - *i£$xr , ("»> 

and for large L 4 but any t 

(0(t)) = (0|O|0) + O (e-OEi-So)^ . (114) 
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One has |/) = e _l< ^ |i) if the fields are periodic up to phases, rotations in flavor, 
and discrete symmetry operations. For example, to set |/) = with no non- 
trivial phase, one chooses periodic boundary conditions for bosonic variables 
and emiz-periodic for fermionic variables: 

A^Li, x) = +A fl (0, x), ^{Li, x) = -V(0, as). (115) 

The minus sign for fermions follows from Fcrmi-Dirac statistics. The sum over 
initial state \i) is achieved by integrating over the fields at x A = 0. 

Let us focus for a while on the boundary conditions in Eq. (115). Then 
one can re-write Eq. (113) as 

TrOe~^ i4 
Tre _ffi4 

Thus, the Euclidean functional integrals give thermal expectation values for a 
quantum system with temperature 

T=(k B L 4 )-\ (117) 

where ks is Boltzmann's constant. Calculations with non-zero temperatures 
T ~ A/fcs are of considerable interest in nuclear physics and astrophysics. 
When the temperature (and also the chemical potential for baryons) is in this 
range, one expects a phase transition. The hot, dense phase is called the 
quark-gluon plasma, and it is thought to exist in heavy-ion collisions 167 and 
in quark stars, 168 which are similar to neutron stars, but with large enough 
gravitational attraction to crush the neutrons into the quark-gluon plasma. 

Lattice calculations have been used to compute the critical temperature 
of the phase transition. In the pure gauge theory, this was one of the first 
large-scale calculations without guidance from experiment. 169 ' 170 The critical 
temperature for QCD with 2 and 4 flavors has also been computed. 171-173 Over 
the last several years, many issues in QCD thermodynamics have been studied, 
such as the equation of state of the quark-gluon plasma. 174 For further reading 
on this rich application of lattice QCD, there are several recent reviews. 175 

For applications of lattice QCD motivated by particle physics, we are in- 
terested in zero temperature. In practice, L4 is large but not infinite, so T is 
small but non-zero. From Eq. (114), one sees that the corrections are small: 
the energy gap between a state and the vacuum is nothing but the energy 
of a particle, as conventionally defined. In quenched calculations, E\ — Eq is 
the glueball mass. 9 With dynamical light quarks, the energy gap will be the 

9 One can consider quenched QCD to contain qq mesons (and qqq baryons) in the Hilbert 
space. But then one must also include qq mesons, etc. The contributions of genuine hadrons 
and the ghosts cancel in Eq. (113). 
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"pion" mass, that is, the mass of the lightest pseudoscalar mps- The expo- 
nential suppression means that to reduce the thermal correction below 0.1% 
one needs m PS L^ > 7 or L 4 > 4 fm for m PS ~ 350 MeV. 

The discussion of Eqs. (108)-(116) assumed continuous (Euclidean) time. 
For discrete time, the analysis stays the same, but the relation between the 
Lagrangian and the Hamiltonian is a little different. The essential idea is as 
follows. There is, perhaps not surprisingly, a natural time-evolution operator, 
called the transfer operator (or, more often, the transfer matrix). 1 ' 10 - 164 " 166 It 
propagates states in Hilbert space forward one unit of Euclidean time: 

\n{ Xi + a))=T\n{ Xi )), (118) 

where |fi) is any state in the Hilbert space. For details on the relationship be- 
tween the lattice action and the transfer matrix, and how to define the Hilbert 
space, the reader should consult one of the textbooks. 12-14 The Hamiltonian H 
is then defined by solving 

T = exp(-aiJ). (119) 

Indeed, all the masses considered in Sec. 2 are eigenvalues of this Hamiltonian. 
The maximal eigenvalue To is simply e~ aE °, where Eq is interpreted as the 
vacuum energy. Equation (119) provides an adequate and useful definition 
of H, if T is positive and has several eigenvalues close to the vacuum. Then 
— a -1 ln(T„/To) are just the particle masses and energies. The transfer matrix 
is indeed positive for simple actions, like the Wilson action. For improved 
actions there can be non-positive eigenvalues; unless the action is pathological 
(and, thus, not an improvement) such eigenstates have energies near the cutoff 
and are, therefore, not physical. 

7.2 Finite spatial volume 

This subsection discusses corrections from the finite spatial volume. There are 
two classes of effects: polarization effects and scattering effects. The former 
arise because interacting particles are surrounded by a cloud of virtual parti- 
cles. In a periodic volume, these virtual particles can propagate "around the 
world." Except when the cloud contains relatively light particles with Lm < 1, 
these are exponentially suppressed. The scattering effects are much more in- 
teresting. The boundary influences the energies of two-particle states in a way 
that is connected to the phase shifts. As a consequence, calculations of the 
L dependence of these energies yields valuable information, such as scattering 
lengths, resonance widths, and so on. We shall aim for a somewhat simplified 
explanation. A more mathematical review is in van Baal's chapter. 176 

We consider the size L of the box to be large enough that the finite volume 
does not alter the structure of the hadrons, namely LA 3> 1. Then finite- 
volume effects are at long distances (by definition) , so the appropriate degrees 
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of freedom, in QCD, are hadrons. As long as the "pion" is not too light, namely 
Lmps 3> 1, the correct effective quantum field theory includes a massive field 
for each hadron. Unlike the effective theories of Sees. 4-6, there is no split 
here into a leading term and small corrections. Instead, the hadronic field 
theory provides a general parametrization that is consistent with analyticity, 
unitarity, and symmetries. The utility of such a theory for describing finite- 
volume effects was first pointed out by Luscher, 177 ' 178 and most of the results 
sketched below are due to him. 177-180 

For the polarization effects, Luscher used an all-orders skeleton expansion 
for self-energies and vertex functions. The main features can be appreciated at 
the one-loop level. For simplicity, let us consider only two particles, a "nucleon" 
and a "pion" , with effective Lagrangian 

C = -N (0 + m QN ) N-\ {d^f - \mlir 2 + ynNj 5 N. (120) 

The Yukawa interaction leads to the nucleon self-energy 

y (n \ = .p f dk i 1 V i{j>+p)+m m 
MW V J 2tt (k*+m*)[(p + k)*+m* N ]' [ ' 

where the spatial momentum is discrete in a finite box. For periodic boundary 
conditions, tt(x + Lei) — n(x), components of the spatial momentum satisfy 

e ik iL =1=>k . = _ Vu (122) 

where Vi is an integer. Our aim is to exhibit the difference between the mode 
sum in Eq. (121) and an integral over all k. 

The first step is to use an exponential representation of the propagators. 
Then the self energy becomes 

roc rl p 11 -i 

Sl(p) = V 2 J a d P I J -^tfT, W +® + 
Now the three sums over spatial momenta are disentangled to the form 



Vi— oc 



Ar,s) = \ ]T e-^+ 2s ^, (124) 

v% — — OO 

where, in our case, r = p(2tt/L) 2 and s — Si = xpiL/2u for each i. The terms 
with kiji in Eq. (123) require 



U—OO 



2r ds 



(125) 



•50 



Handbook of QCD / Volume 4 



The sum J can be rearranged with the help of the Poisson resummation formula 

1 e px 2 p 2 t +2^ f-n 2 L 2 \ 
I( p (2tt/L) ,Si) = —= 2^ cos(nxpiL) exp f — j. 

V yr n __ GO \ A / 



(126) 



When L is large, the term with n = dominates; it is exactly the same as the 
the integral over all fc,. The leading finite- volume correction comes from the 
terms with n = ±1. Thus, 



1 P px 2 Pi / 

T(p(27r/£) 2 ,s 4 ) = p= F (l + 2cos(xp lJ L)e- i / 4 " (127) 

V47T ^JP \ ' 

= J g e -P(fe 2 +2^fe) (i + 2cos(x Pj i)e- L2 /4^ (128) 

up to terms suppressed by e~ L l p . 

To derive the finite- volume correction to the nucleon mass, it is enough to 
set p on shell: p4 = 2m 7v74, p = 0. Then, 



P OO P 1 

Jo Jo 



2 1 
(47r)<V 2 

x (l + 2e- L2 / 4p y~\ (129) 



in d dimensional spacetime. Eq. (129) is remarkable for several reasons. It 
naturally lends itself to dimensional regularization of ultraviolet divergences, 
which, when d > 4, arise from the lower limit of the integration over p. But 
in the finite- volume corrections the exponential factors e~™ L / 4p suppress this 
region, faster than any power of p. Therefore, renormalization of ultraviolet 
divergences is carried out the same way as in infinite volume. In Eq. (129), the 
divergence is removed by mass renormalization. The asymptotic finite-volume 
correction can be evaluated using the method of steepest descent to carry out 
the integration over x and p. Assuming m 2 < 2m 2 N , 

m N (L) - mjv(oo) = - S L = -■^--^f^— r~ (130) 

87T m z N L 

where p? = m 2 (l — m 2 ./ 4:171%). The finite volume mass shift is less then 0.1% 
if Lp > 7. 

The L dependence holds for any stable particle and to all orders in the 
interaction. 177 (Resonances, such as the p meson and the A baryon, are another 
matter; see below.) If both particles in the self energy have the same mass m, 
then p = \^3m/2. In the case of the physical pion and nucleon, m n <C 2niN, 



Lattice QCD 



51 



so /j, w Tot,.. Even if the lightest pseudoscalar has mps = 400 MeV, which 
corresponds to a light quark mass m s /3, then /i w mpg. 

For matrix elements with at most one hadron in the final state, one can 
carry out a similar analysis. The basic structure remains the same: use an 
exponential representation of the propagators and then evaluate the mode 
sums with the Poisson formula. In this way, it follows that the finite volume 
corrections are exponentially suppressed (as long as Lm 3> 1 for all particles) . 

Equation (130) holds only if Lm^ is large. When Lm-n < 1, the integration 
over x in Eq. (129) is qualitatively different. For m n — 0, the x integration is 
elementary, leading to 

3y 2 1 

m N {L) - TOjv(oo) = L 2 mN (for 771^ = 0). (131) 

Thus, as one would expect, finite- volume corrections are larger when there is a 
cloud of massless particles. Equation (131) is a simple example of finite- volume 
chiral perturbation theory, which can be used rather generally for both finite 
volume and non-zero temperature. 181 " 183 Section 6 noted that the logarithms 
of chiral loops are difficult to demonstrate in numerical data. That is unfor- 
tunate, because chiral behavior is a good a posteriori check of algorithms for 
full QCD; indeed, the check may be easier to understand than the algorithms 
themselves. It may, therefore, prove worthwhile to hold L artificially small, 
and test for finite-volume Goldstone-boson effects. 

Before turning to finite-volume effects in scattering states, let us consider 
how to extract hadron masses and matrix elements with physical light quark 
masses and infinite volume, from numerical data with practical light quark 
masses and finite volume. In practice, the lightest pseudoscalar is indeed mas- 
sive, because of the chiral slowing down of the numerical algorithms, Eq. (21). 
For the sake of illustration, let us assume the lightest pseudoscalar has a mass 
tops > 400 MeV. A tractable situation has L = 2 fm, L4 = 2L. The expo- 
nential factor e ~ mpsL then falls below 2% (and e - mpsL * below 0.05%). Up 
and down quarks are reached by fitting data with to s /3 < m q < m s to the 
guide given by infinite-volume chiral perturbation theory. Although the finite- 
volume corrections from Eq. (130) are present in the data, they induce only a 
small bias in the fit. The fit parameters inherit the bias, but when reconstitut- 
ing the desired hadronic property from the fit parameters and physical quark 
masses, the finite- volume effect remains under control. 

Finite-volume effects are most interesting for scattering states, including 
resonances, such as the p meson and the A nucleon. Using models Wiese 184 and 
DeGrand 185 showed that it would be difficult to determine resonance masses 
in a finite volume. The main difficulty is that a resonance has the same quan- 
tum numbers as multi-particle states. Intuitively, the latter depend strongly 
on L, because the natural unit of relative momenta is 2n/L, but a resonance's 
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mass should be nearly independent of L. Hence, there are many would-be 
level crossings. Of course, levels cannot cross as a parameter (here L) varies 
smoothly; instead the states interact and the levels push each other apart. It 
is plausible that this phenomenon is sensitive to details of the dynamics, but 
details are not immediate. 

A genuine breakthrough came in Liischer's 1991 paper, 179 which derived 
a non-perturbative relation between two-particle energies and phase shifts. In 
a finite volume, the energy spectrum is discrete. The main insight is that 
the phase shift is generated at relatively short range, whereas the boundary 
conditions that render a discrete spectrum are imposed at very long distances. 
Therefore, the two aspects of the problem can be attacked separately. Fur- 
thermore, it is clear from this consideration that it is the infinite volume phase 
shift that should arise. 

The relation between the spectrum and phase shifts has been proven for 
elastic scattering. In the inelastic region, conceptual problems with the defini- 
tion of finite-volume n-particle states (n > 2) have not been worked out. For 
pion-pion scattering, the elastic region is 2m^ < E < Am n ; for pion-nucleon 
scattering, mN+m v < E < mjv + 3m T . Particularly in numerical lattice QCD, 
where the pion is artificially heavy, these limits apply to the p and A. Indeed, 
in many examples we suppose that mps > 300 GcV; for p and A resonance 
properties, one should reduce the pseudoscalar mass further. 

In infinite volume, the scattering phase shift arises in connecting the short- 
and long-range parts of the two-body wave function. In QCD long range means 
distances greater than tq ~ 1-2 fm. Thus, the phase shift encodes the part 
of the confining interactions of quarks and gluons that is "remembered" after 
hadrons have scattered. At very large separation, the only boundary condition 
on the two-body wave function is that its scattered component fall off like a 
spherical wave. All energies are possible and chosen by the initial conditions. 
Thus the scattering amplitude (which depends on the phase shifts) is the most 
interesting object, and experiments measure cross sections. 

In a periodic finite volume, several features are radically different. Recall 
that one should take the time extent L 4 is at least as big as L, to set the 
temperature to zero. In that case, the two "scattering" particles either never 
separate very far (if their relative momentum is small) or encounter each other 
more than once (if large). So, a straightforward scattering problem does not 
arise; scattering amplitudes and cross sections are not natural observables. On 
the other hand, the two-body wave function has to be periodic, 

Mr + e t L) =Vfe(r), (132) 

where r is the relative separation of the two particles, and the label k corre- 
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sponds to the relative momentum. The energy of the two-body state is 



E = 2^m 2 + k 2 , 



(133) 



when both particles are assumed to have the same mass to. For fixed L, 
Eq. (132) allows only certain values of k. Because the boundary condition is 
a long-distance effect, the only aspect of the interaction that can enter is the 
set of phase shifts. The allowed values can be compactly (but impenetrably) 
summarized as the solutions to 179 



det 



U(kL/2n) = 0, 



(134) 



where the matrix indices are azimuthal and magnetic angular momentum quan- 
tum numbers I and to, f/j m ,j'm'- The phase shift matrix is 



= 2j<5(fc) 



ImJ' 



(135) 



where Si(k) is the phase shift in the i-th partial wave, and (in an unfortunate 
clash of notation) Sw is the Kronecker delta. 

The most important feature of Eq. (134) is that U(q) is built from com- 
pletely kinematical, L independent functions. 179 Moreover, U splits into blocks 
along the diagonal, with a block for each irreducible representation of the cubic 
group. Finally, if, in each block, higher partial waves have small phase shifts, 
then Eq. (134) reduces to a simple equation. 

To illustrate this structure, let us consider the case of the p meson, con- 
sidering the practical situation m p < Am^ . First let us neglect all phase shifts 
with / > 4. Then Eq. (134) becomes 



e 2«i(fc)_ Un 



-«13 

e 2 *»< fc ) -U33 



0. 



(136) 



This is tractable, because the information needed to get Uij has been 
tabulated. 179 It is even simpler if one can neglect <5 3 . 180 Then 



<5i(fc) = nir - 4>{kL/2ir), 



(137) 



where 



tan (j)(q) = — 
Zoo 



7r 3 /y 

£oo(l;<7 2 ) 
1 

47T 



(138) 
(139) 
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with the sum over three-vectors of integers. The zeta function Z o (and sev- 
eral cousins Zi m ) arc defined by the sum for Re 2s > I + 3 and by analytic 
continuation for s = 1. <p(q) is a monotonically increasing function of q; as one 
pieces together the phase shift, the integer n is chosen so that Si(k) is smooth. 

To see how this works, let us outline how to obtain the p mass and width. 180 
The first step is to compute, via Monte Carlo, the 1=1 two-pion ener- 
gies E, as a function of L. Via Eq. (133) each energy yields a momentum 
k = \ \JE 2 — Am 2 , and then the phase shift Si(k) is obtained from Eq. (137). 
The p resonance appears at momentum k p , such that 5\{k p ) = ir/2. Then, by 
studying the L (and hence k dependence) around the resonance, the effective 
range formula 

k 3 Akl{k 2 n -k 2 ) 

¥ cot f !(<»)= p[p ' (140) 

can be used to determine the width. The resulting {m p ,T p ) would correspond 
to an unphysical case, with too large quark masses. From chiral perturbation 
theory, we expect m p to depend mildly on quark masses. Similarly, if the decay 
is modeled with the interaction 

C P1T7T = g P ^s abc p a pTr b d t "TT c , (141) 



then xPT suggests that g plxlx depends only mildly on the mass. The width 

n 2 k 3 

T p = (142) 



derived from C pm r is 



p 



with k 3 accounting for the broadening of the resonance as phase space opens 
up. Thus, one could determine the physical p meson width by extrapolating 
g pn7r to physical quark mass, and then using Eq. (142). 

The relationship between the energy levels of scattering states and phase 
shifts extends to weak non-leptonic decays, again in the elastic region. 186 Thus, 
it is not a conceptual problem to calculate the phase shifts in K decays. Un- 
fortunately, for for inelastic decays, such as B meson decays, it remains an 
unsolved problem. 

Often an even bleaker picture is painted, based on a superficial under- 
standing a theorem of Maiani and Testa. 187 The theorem assumes an infinite 
volume and is, thus, relevant only to extremely large volumes, where there the 
level spacing is so small that it is not possible, because of statistical errors, to 
resolve them. Then one simply observes that Euclidean correlation functions 
are real, so phases cannot appear without analytic continuation to Minkowski 
space. Although possible in principle, 9 this is again not practical when the 
correlation functions are computed only at discrete values of time separation. 
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As with many theorems, the way out is to relax one of the hypotheses: in this 
case, the assumption of infinite volume. 

Although the study of finite volume effects, in the elastic region, is on 
a sound conceptual footing, there are several technical challenges. First, as 
mentioned above, the quark masses must be somewhat smaller than usually 
used. Second, many independent Monte Carlo runs, with varying L but other 
parameters fixed, are needed. Third, the quenched approximation is not suited 
to these calculations at all, because the would-be rf produces unphysical terms 
in the 1/L expansion. 188 Finally, small statistical errors are needed, to trace 
out the level-crossing regions, especially when the resonance is narrow. Further 
practical aspects have been considered by Lin et al. 189 

8 Top Ten Trends 

When the Editor invited me to contribute a chapter for the Handbook of QCD, 
one possibility that was raised was a list of the top ten results from lattice QCD. 
One version of this idea, which was quickly discarded, would have been a list 
of specific numerical results. Such lists are quickly out of date, and are better 
obtained from recent from lattice conferences. 

Nevertheless, it is tempting to list some of the main trends. In contrast to 
the rest of the chapter, which focused on practical aspects, the following topics 
should enjoy broad interest and influence one's thinking about QCD. 

1. Hadron spectrum: Ab initio calculations of the hadron spectrum were one 
of the early motivations for numerical lattice QCD. The results of a re- 
cent calculation 190 are shown in Fig. 6, comparing quenched and nj = 2 
calculations. The quenched approximation does surprisingly well, agree- 
ing with experiment within several percent. For hadrons with strange 
quark, particularly the </> and the O - , the unquenched calculations agree 
even better. For hadrons with only up and down quarks, for example the 
nucleon, it remains difficult to control finite-size effects. 

2. Glueballs must exist: When QCD was first put forward, it seemed plau- 
sible, but not obvious, that gluons would bind together to form hadrons 
that lie beyond the classification scheme of the quark model. Such states 
are called glueballs and hybrids. Although there is still controversy over 
whether such glueballs have been observed in the laboratory, lattice cal- 
culations give strong evidence that their existence is a prediction of QCD. 
There is even a well-argued case that the /o(1710) is the lowest-lying 
scalar glueball, based on calculations of the masses and coupling to pseu- 
doscalar meson pairs. 191 

3. Confinement explained: An early result of the strong coupling expansion 
of lattice gauge theory was the so-called "area law" for Wilson loops. A 
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Figure 6: The hadron spectrum. 190 Solid (open) point show the results of unquenched 
(quenched) lattice QCD. Blue (red) points tune the strange quark's mass to the K (cj>). 



Wilson loop is the line integral of the gauge potential around a closed 
curve C, or W(C) — P exp (J c A ■ dxj , where P denotes path ordering. 
For a R x T rectangle, a short calculation in the strong coupling limit 
shows 1 

W(R,T)~\^j , (143) 

namely, that lnVF(C) is proportional to the area enclosed by C. For 
large T, bxW(R,T) ~ -V(R)T, where V(R) is the potential energy of 
the gluon field between two static sources of color. Thus, the area law 
shows a linearly rising potential, which confines quarks. It is amazing 
that confinement emerges so simply. 

Numerical calculations confirm that the linear rise persists for non- 
Abelian gauge theories into the weak (bare) coupling regime, i.e., towards 
the continuum limit. (It goes away for the Abelian U(l) lattice gauge the- 
ory.) A recent survey is in the review of Bali. 192 The lattice calculations 
trace out unmistakably a potential that is Coulombic at short distances 
(as expected from asymptotic freedom) and linear at large distances (as 
needed for confinement). 

Insight into confinement is not limited to the heavy-quark potential. 
From a chromoelectric perspective, studies of the gluonic flux between 
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static sources arc starting to paint an even more detailed picture. 193 
From a complementary perspective, several ideas for the condensation of 
chromomagnctic monopoles have been studied. 194 

4. NRQCD and HQET: It is worth recalling that both non-relativistic QCD 
(NRQCD) and heavy-quark effective theory (HQET) were invented, in 
part, to facilitate lattice calculations of heavy-quark systems. 99,101 With 
continuum ultraviolet regulators, instead of a lattice, these effective field 
theories have become powerful tools for studying, for example, quarko- 
nium production in high-energy collisions, or semi-leptonic decays of B 
mesons. They are important for understanding the CKM matrix, and 
QCD itself. 

5. The light quarks are very light: There is now considerable evidence from 
lattice QCD that the strange quark has a mass lower than previously 
thought. The up and down quarks arc correspondingly lighter, as ex- 
pected from chiral symmetry: (m u + m<j)/2m s ~ m 2 /2m 2 K rts 1/25. A 
recent calculation with dynamical quark loops finds the MS mass 190 

m s (2 GeV) = 88± 4 MeV. (144) 

This was not an easy conclusion to reach, because lattice spacing effects 
and quenching effects are both large. If either is not taken into account, 
a larger estimate of fh s is obtained. One implication of this result is 
that very little of the mass of everyday objects comes from the Higgs 
mechanism of the electroweak theory. Almost all of it comes from gluons. 

6. QCD has a phase transition: There are theoretical grounds to expect 
that pure gauge theory has a deconfinement transition, and that in QCD 
the high-temperature phase restores chiral symmetry. (See Part 8 of 
Vol. 3 of the Handbook.) But without numerical lattice QCD it is im- 
possible to compute the critical temperature of the phase transition(s). 
The calculation of the pure gauge transition temperature was one of the 
first large-scale numerical simulations to yield, persuasively, a physically 
interesting result that was not known from experiment. 

7. a s from hadrons agrees with a s from high-energy scattering: One of the 
entries for a s in the Review of Particle Properties 195 comes from lattice 
calculations of charmonium and bottomonium spectra. It lies within the 
range of the other determinations, which come from high-energy scatter- 
ing. The agreement shows that the QCD of confinement is also the QCD 
of partons: QCD works at long and short distances. 

8. Determination of the CKM matrix: Lattice calculations of matrix ele- 
ments in B, D, and K decays are now seen as essential to testing the 
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standard CKM picture of flavor and CP violation. 96 The most appealing 
aspect of the CKM mechanism is that it is very predictive, relating many 
CP conserving and CP violating processes back to four parameters, only 
one of which violates CP. With measurements of CP conserving pro- 
cesses, such as semi-leptonic decays and neutral meson mixing, and cal- 
culations of the corresponding form factors and mixing amplitudes, one 
can predict the strength of CP violation. Once unquenched QCD calcu- 
lations of these quantities mature, this area could well become the most 
important contribution of lattice QCD. 

9. Lattice chiral symmetry: The last few years have witnessed remarkable 
progress of understanding chiral symmetry in lattice gauge theory. 79-81 
In the context of QCD, the contribution is significant: a tool for more rig- 
orous derivations of well-known soft-pion theorems, 196 or an ingredient in 
calculations of kaon matrix elements relevant to CKM phenomenology. 86 
The contribution of this development of lattice gauge theory to elec- 
troweak physics and to extensions of the Standard Model promises to be 
much greater. 

10. Perturbation theory is universal: A somewhat underappreciated, but nev- 
ertheless important, achievement is the proof that perturbation theory 
is universal. 68 After renormalization, the a — > limit of perturbative 
QCD is the same for all lattice actions. Since lattice field theory main- 
tains rigorous control over the ultraviolet at all steps of the analysis, it 
puts more familiar (i.e., continuum) ultraviolet regulators on a sounder 
footing, because they are all connected to lattice gauge theory through 
BPHZ renormalization. 

9 Summary: Computational vs. Theoretical Physics 

Lattice gauge theory is a broad subject, with many applications beyond QCD, 
especially now that chiral fermions are better understood. The lattice pro- 
vides a well-defined mathematical foundation for quantum field theory, which 
remains one of its strongest attractions. In particular, it permits, in least in 
principle, real non-perturbative calculations. 

Especially in the realm of QCD, one technique for calculation is most 
prominent, namely computing the functional integral numerically, via Monte 
Carlo methods. Because we live in 3+1 dimensional spacetime, and because 
real hadrons have several scales (m q , A, toq), the problem is far from easy. 
When discussing numerical lattice QCD, usually only the computational as- 
pects that are emphasized, perhaps overemphasized. 

Computing is necessary, but not sufficient. In addition to the physical 
scales of hadrons, the numerical technique has an ultraviolet cutoff (from the 
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lattice) and an infrared cutoff (from a finite spacetime volume). A whole branch 
of theoretical physics has been developed to understand the cutoff effects, and 
how to get the most out of scarce computer resources. Similarly, for practical 
reasons one must cope with light quark masses that are not so light, and 
heavy quark masses near the ultraviolet cutoff. Here tools are available from 
continuum QCD, but they require some redesign to be usefully deployed. 

Even if the numerical solution of QCD were accessible solely by brute force, 
it would be wasteful to rely on brute force alone. An example, discussed at the 
end of Sec. 4, is the control of cutoff effects. It is much more efficient to run 
at several lattice spacings than to put all computer resources onto the finest 
conceivable lattice. A similar comment applies to the light quark mass, but 
here, at least, it is completely routine to do so. 

In the past decade computing has been essential to test several new ideas, 
such as faster algorithms, methods for heavy quarks, and non-perturbative 
rcnormalization. Most of the progress in lattice QCD has come, however, from 
the ideas and not from Moore's Law, which observes that computer power 
per dollar doubles every 18 months. Because we live in four dimensions, and 
because the algorithms slow down, one can double the lattice spacing (at fixed 
computer cost) only every 6-8 years. Increases in computing are essential, but 
not sufficient. 

It is likely that this trend will continue. For example, non-perturbative 
renormalization is helpful for light quarks, but remains immature for heavy 
quarks. For kaon physics relatively new, chirally symmetric, methods for the 
quarks are promising. But they are even more computationally intensive. One 
of the methods, domain wall fermions, adds a fifth dimension to the lattice 
and to the scaling laws in Sec. 3. The other methods are similarly difficult. 
Thus, especially here, the need for ideas to reduce the computational burden 
(for fixed error bar) is great. 

Even if no new ideas emerge plenty of traditional theoretical physics is 
needed to understand numerical lattice gauge theory and, therefore, QCD. Be- 
cause effective field theories provide tools to reduce the uncertainties, there 
is always more work to do: calculating the next order in chiral perturbation 
theory, or programming the next set of operators to reduce discretization ef- 
fects. This work is necessary, especially when the results of lattice QCD are 
needed to understand flavor or collider physics. One should assume, however, 
the there will be new ideas, and that the theoretical and computational sides 
will continue to inform each other. 

Acknowledgments 

I would like to thank Misha Shifman for the invitation to write this chapter, 
and for his patience as it was being written. Over the years, many people 



GO 



Handbook of QCD / Volume 4 



have influenced my views on lattice gauge theory, particularly Peter Lepage, 
Martin Luschcr, and Kenneth Wilson. My collaborators in developing a theory 
of heavy-quark discretization effects, Aida El-Khadra, Shoji Hashimoto, Paul 
Mackenzie, and Tetsuya Onogi, have had a significant influence on my way of 
thinking. While writing this chapter, I benefitted from a conversation on chiral 
perturbation theory with Maarten Golterman and Stephen Sharpe. Finally, I 
would like to thank the Aspen Center for Physics and Zoltan Ligeti for provid- 
ing an opportunity to rehearse the material before a live (and lively!) audience. 
Fcrmilab is operated by Universities Research Association Inc., under contract 
with the U.S. Department of Energy. 

1. K. G. Wilson, Phys. Rev. D10, 2445 (1974). 

2. J. D. Bjorkcn and E. A. Paschos, Phys. Rev. 185, 1975 (1969). 

3. H. Fritzsch, M. Gell-Mann, and H. Lcutwylcr, Phys. Lett. B47, 365 
(1973). 

4. D. J. Gross and F. Wilczek, Phys. Rev. Lett. 30, 1343 (1973). 

5. H. D. Politzer, Phys. Rev. Lett. 30, 1346 (1973). 

6. J. J. Aubert et al, Phys. Rev. Lett. 33, 1404 (1974). 

7. J. E. Augustin et al, Phys. Rev. Lett. 33, 1406 (1974). 

8. F. J. Wegner, J. Math. Phys. 12, 2259 (1971). 

9. J. Glimm and A. Jaffc, Quantum Physics: A Functional Integral Point 
of View, 1st. edition (Springer, New York, 1981); 2nd. edition (Springer, 
New York, 1987). 

10. K. G. Wilson, in New Phenomena in Subnuclear Physics, edited by A. 
Zichichi (Plenum, New York, 1977). 

11. K. G. Wilson and J. B. Kogut, Phys. Rept. 12, 75 (1974). 

12. M. Creutz, Quarks, Gluons, and Lattices (Cambridge University, Cam- 
bridge, UK, 1983). 

13. H. J. Rothe, Lattice Gauge Theories: An Lntroduction, 1st. edition 
(World Scientific, Singapore, 1992); 2nd. edition (World Scientific, Sin- 
gapore, 1997). 

14. I. Montvay and G. Minister, Quantum Fields on a Lattice (Cambridge 
University, Cambridge, UK, 1994). 

15. J. B. Kogut, Rev. Mod. Phys. 55, 775 (1983). 

16. A. Hasenfratz and P. Hasenfratz, Annu. Rev. Nucl. Part. Sci. 35, 559 
(1985). 

17. A. S. Kronfeld and P. B. Mackenzie, Annu. Rev. Nucl. Part. Sci. 43, 793 
(1993) [hep-ph/9303305]. 

18. S. Aoki, in Proceedings of the XIX International Symposium on Lepton 
and Photon Interactions at High Energy, edited by J. A. Jaros and M. E. 
Peskin, eConfC 990809, 657 (2000) [hep-ph/9912288]. 

19. R. D. Kenway, in ICHEP 2000: Proceedings, edited by C. S. Lim and 



Lattice QCD 



61 



T. Yamanaka (World Scientific, Singapore, 2001) [hep-ph/0010219]. 

20. From Actions to Answers, edited by T. DeGrand and D. Toussaint 
(World Scientific, Singapore, 1990). 

21. M. Liischcr, in Fields, Strings, and Critical Phenomena, edited by E. 
Brezin and J. Zinn-Justin (North- Holland, Amsterdam, 1990). 

22. A. S. Kronfeld, in Perspectives in the Standard Model, edited by R. K. 
Ellis, C. T. Hill, and J. D. Lykken (World Scientific, Singapore, 1992). 

23. S. R. Sharpe, in CP Violation and the Limits of the Standard Model, 
edited by J F. Donoghue (World Scientific, Singapore, 1995) [hep- 
ph/9412243]. 

24. R. Gupta, in Probing the Standard Model of Particle Interactions, edited 
by R. Gupta, A. Morel, E. De Rafael, F. David (Elsevier, Amsterdam, 
1999) [hep-lat/9807028]. 

25. M. Liischcr, in Probing the Standard Model of Particle Interactions, 
edited by R. Gupta, A. Morel, E. De Rafael, F. David (Elsevier, Ams- 
terdam, 1999) [hep-lat/9802029]. 

26. G. Minister and M. Walzl, in Phenomenology of Gauge Interactions, 
edited by D. Graudenz and V. Markushin (URL http:/ /ltpth.web.psi.ch/ 
/zuoz/zuoz2000/zuoz2000proc.htm) [hep-lat/0012005]. 

27. T. DeGrand, in Flavor Physics for the Millennium, edited by J. L. Rosner 
(World Scientific, Singapore, 2001) [hep-ph/0008234]. 

28. W. T. Giele, S. A. Keller, and D. A. Kosower, hcp-ph/0104052. 

29. M. Gockeler et al, Phys. Rev. D53, 2317 (1996) [hep-lat/9508004]; 

30. M. Guagnclli, K. Janscn, and R. Pctronzio, Phys. Lett. B493, 77 (2000) 
[hep-lat/0009006]; 

31. D. Dolgov et al. [LHPC and SESAM Collaborations], hcp-lat/0201021. 

32. K. Janscn, in ICHEP 2000: Proceedings, edited by C. S. Lim and T. Ya- 
manaka (World Scientific, Singapore, 2001) [hep-lat/0010038]. 

33. M. Di Pierro, "From Monte Carlo Integration to Lattice Quantum Chro- 
modynamics: An Introduction," hep-lat/0009001 [FERMILAB-FN-699] . 

34. H. Neuberger, Phys. Lett. B417, 141 (1998) [hcp-lat/9707022]; B427, 
353 (1998) [hep-lat/9801031]. 

35. L. Susskind, Phys. Rev. D16, 3031 (1977); T. Banks, J. B. Kogut, and 
L. Susskind, ibid. D13, 1043 (1976). 

36. S. Aoki et al. [JLQCD Collaboration], Phys. Rev. D65, 094507 (2002) 
[hcp-lat/0112051]. 

37. S. Duanc, Nucl. Phys. B257, 652 (1985). 

38. A. Ukawa and M. Fukugita, Phys. Rev. Lett. 55, 1854 (1985). 

39. G. G. Batrouni et al, Phys. Rev. D32, 2736 (1985). 

40. J. C. Sexton, private communication; A. C. Irving and J. C. Sexton, 
Phys. Rev. D55, 5456 (1997) [hep-lat/9608145]. 

41. E. Marinari, G. Parisi, and C. Rebbi, Nucl. Phys. B190, 734 (1981). 



62 



Handbook of QCD / Volume 4 



42. D. Weingarten, Phys. Lett. B109, 57 (1982). 

43. J. Sexton and D. Weingarten, Phys. Rev. D55, 4025 (1997). 

44. A. X. El-Khadra, G. Hockney, A. S. Kronfeld, and P. B. Mackenzie, Phys. 
Rev. Lett. 69, 729 (1992). 

45. P. B. Mackenzie, Nucl. Phys. B Proc. Suppl. 34, 400 (1994). 

46. C. T. H. Davies et al, Phys. Rev. Lett. 73, 2654 (1994) [hep-lat/9404012]. 

47. B. J. Gough et al, Phys. Rev. Lett. 79, 1622 (1997) [hep-ph/9610223]. 

48. S. Hashimoto et al, hcp-ph/0110253; J. N. Simone et al, Nucl. Phys. B 
Proc. Suppl. 83, 334 (2000) [hep-lat/9910026]. 

49. A. Ali Khan et al. [CP-PACS Collaboration], Phys. Rev. D64, 034505 
(2001) [hep-lat/0010009]. 

50. A. Ali Khan et al. [CP-PACS Collaboration], Phys. Rev. D64, 054504 
(2001) [hcp-lat/0103020]. 

51. C. Bernard et al. [MILC Collaboration], Nucl. Phys. B Proc. Suppl. 83, 
289 (2000) [hep-lat/9909121]; 94, 346 (2001) [hep-lat/0011029]. 

52. V. Gimcncz, L. Giusti, G. Martinelli, and F. Rapuano, JHEP 0003, 018 
(2000) [hep-lat/0002007]. 

53. S. Coleman and R. E. Norton, Nuovo Cim. 38, 438 (1965). 

54. A. S. Kronfeld, Prog. Theor. Phys. Suppl. Ill, 293 (1993) [hep-lat/ 
9205008]. 

55. S. Gottlieb, W. Liu, D. Toussaint, R. L. Renkcn, and R. L. Sugar, Phys. 
Rev. D35, 2531 (1987). 

56. S. Duane, A. D. Kennedy, B. J. Pendleton, and D. Roweth, Phys. Lett. 
B195, 216 (1987). 

57. M. Liischcr, Nucl. Phys. B418, 637 (1994) [hep-lat/9311007]. 

58. G. P. Lepage et al, Nucl. Phys. B Proc. Suppl. 106, 12 (2002) [hep- 
lat/0110175]. 

59. K. Symanzik, in Recent Developments in Gauge Theories, edited by G. 
't Hooft et al. (Plenum, New York, 1980). 

60. K. Symanzik, in Mathematical Problems in Theoretical Physics, edited 
by R. Schradcr et al. (Springer, New York, 1982); Nucl. Phys. B226, 
187, 205 (1983). 

61. B. Sheikholeslami and R. Wohlert, Nucl. Phys. B259, 572 (1985). 

62. K. Symanzik, Commun. Math. Phys. 18, 227 (1970); 23, 49 (1971). 

63. C. G. Callan, Phys. Rev. D2, 1541 (1970). 

64. M. Liischer and P. Weisz, Commun. Math. Phys. 97, 59 (1985); (E) 98, 
433 (1985). 

65. M. Liischer, S. Sint, R. Sommcr, and P. Weisz, Nucl. Phys. B478, 365 
(1996) [hep-lat/9605038]. 

66. T. Reisz, Commun. Math. Phys. 116, 81 (1988). 

67. T. Reisz, Commun. Math. Phys. 117, 79 (1988). 

68. T. Reisz, Nucl. Phys. B318, 417 (1989). 



Lattice QCD 63 

69. R. Wohlcrt, DESY 87-069 (unpublished, July 1987). 

70. M. Liischer and P. Weisz, Nucl. Phys. B479, 429 (1996). 

71. S. Sint and P. Weisz, Nucl. Phys. B502, 251 (1997). 

72. M. Bochicchio, L. Maiani, G. Martinelli, G. Rossi, and M. Testa, Nucl. 
Phys. B262, 331 (1985). 

73. K. Janscn et al, Phys. Lett. B372, 275 (1996) [hep-lat/9512009]; 
M. Liischer, S. Sint, R. Sommer, P. Weisz, and U. Wolff, Nucl. Phys. 
B491, 323 (1997) [hep-lat/9609035]. 

74. M. Liischer, S. Sint, R. Sommer, and H. Wittig, Nucl. Phys. B491, 344 

(1997) [hep-lat/9611015]. 

75. T. Bhattacharya, R. Gupta, W. Lee, and S. Sharpe, Phys. Rev. D63, 
074505 (2001) [hep-lat/0009038]. 

76. S. Collins, C. T. Davics, G. P. Lepage, and J. Shigemitsu, hep- 
lat/0110159. 

77. S. R. Sharpe, Nucl. Phys. B Proc. Suppl. 34, 403 (1994) [hcp-lat/ 
9312009]. 

78. Y. Luo, Phys. Rev. D55, 353 (1997) [hep-lat/9604025]. 

79. F. Niedermayer, Nucl. Phys. B Proc. Suppl. 73, 105 (1999) [hep-lat/ 
9810026]. 

80. H. Neuberger, Nucl. Phys. B Proc. Suppl. 83, 67 (2000) [hep-lat/ 
9909042]. 

81. M. Liischer, Nucl. Phys. B Proc. Suppl. 83, 34 (2000) [hcp-lat/9909150]. 

82. Y. Shamir, Nucl. Phys. B406, 90 (1993) [hep-lat/9303005]. 

83. U. J. Wiese, Phys. Lett. B315, 417 (1993) [hep-lat/9306003]. 

84. P. Hasenfratz, V. Lalicna, and F. Niedermayer, Phys. Lett. B427, 125 

(1998) [hcp-lat/9801021]. 

85. D. B. Kaplan, Phys. Lett. B288, 342 (1992) [hep-lat/9206013]. 

86. T. Blum and A. Soni, Phys. Rev. D56, 174 (1997) [hep-lat/9611030]. 

87. A. Ali Khan et al. [CP-PACS Collaboration], Phys. Rev. D64, 114506 
(2001) [hep-lat/0105020]. 

88. T. Blum et al. [RBC Collaboration], hep-lat/0110075. 

89. G. Kilcup, R. Gupta, and S. R. Sharpe, Phys. Rev. D57 (1998) 1654 
[hcp-lat/9707006]. 

90. S. Aoki et al. [JLQCD Collaboration], Phys. Rev. Lett. 80 (1998) 5271 
[hcp-lat/9710073]. 

91. M. Liischer and P. Weisz, Nucl. Phys. B266, 309 (1986). 

92. M. A. Nobes, H. D. Trotticr, G. P. Lepage, and Q. Mason, Nucl. Phys. 
B Proc. Suppl. 106, 838 (2002) [hep-lat/0110051]. 

93. S. J. Brodsky, G. P. Lepage, and P. B. Mackenzie, Phys. Rev. D28, 228 
(1983). 

94. G. P. Lepage and P. B. Mackenzie, Phys. Rev. D48, 2250 (1993) [hep- 
lat/9209022]. 



G4 



Handbook of QCD / Volume 4 



95. J. Harada, S. Hashimoto, A. S. Kronfcld, and T. Onogi, Phys. Rev. D67, 
014503 (2003) [hep-lat/0208004]. 

96. M. Beneke, "CP violation and the CKM matrix," hcp-lat/0201011. 

97. For reviews of the calculations themselves see recent (and future) reviews 
at lattice conferences, such as S. Hashimoto, Nucl. Phys. B Proc. Suppl. 
83, 3 (2000) [hep-lat/9909136]; C. Bernard, ibid. 94, 159 (2001) [hep- 
lat/0011064]; S. Ryan, ibid. 106, 86 (2002) [hcp-lat/01 11010]. 

98. W. E. Caswell and G. P. Lepage, Phys. Lett. B167, 437 (1986). 

99. G. P. Lepage and B. A. Thacker, Nucl. Phys. B Proc. Suppl. 4, 199 

(1987) ; B. A. Thacker and G. P. Lepage, Phys. Rev. D43, 196 (1991). 

100. G. P. Lepage, L. Magnea, C. Nakhleh, U. Magnea, and K. Hornbostel, 
Phys. Rev. D46, 4052 (1992) [hcp-lat/9205007]. 

101. E. Eichtcn, Nucl. Phys. B Proc. Suppl. 4, 170 (1987). 

102. E. Eichtcn and B. Hill, Phys. Lett. B234, 511 (1990); B240, 193 (1990). 

103. E. Eichtcn and B. Hill, Phys. Lett. B243, 427 (1990). 

104. B. Grinstein, Nucl. Phys. B339, 253 (1990). 

105. H. Gcorgi, Phys. Lett. B240, 447 (1990). 

106. M. E. Luke, Phys. Lett. B252, 447 (1990). 

107. G. T. Bodwin, E. Braatcn, and G. P. Lepage, Phys. Rev. D46, 1914 
(1992) [hcp-lat/9205006]. 

108. A. X. El-Khadra, A. S. Kronfcld, and P. B. Mackenzie, Phys. Rev. D55, 
3933 (1997) [hcp-lat/9604004] . 

109. A. S. Kronfcld, Phys. Rev. D62, 014505 (2000) [hep-lat/0002008]. 

110. J. Harada, S. Hashimoto, K.-I. Ishikawa, A. S. Kronfcld, T. Onogi, and 
N. Yamada, hep-lat/0112044. 

111. J. Harada, S. Hashimoto, A. S. Kronfcld, and T. Onogi, hep-lat/0112045. 

112. A version of this argument is given in S. Aoki, Y. Kuramashi, and S.- 
I. Tominaga, hep-lat/0 107009. This work misstates the origin of the 
\J mi/m,2 in Eq. (68), and it also counts redundant directions incorrectly. 

113. B. P. G. Mertens, A. S. Kronfeld, and A. X. El-Khadra, Phys. Rev. D58, 
034505 (1998) [hep-lat/9712024]. 

114. M. B. Gavela, L. Maiani, S. Petrarca, G. Martinelli, and O. Pene, Phys. 
Lett. 206B, 113 (1988). 

115. C. Bernard, T. Draper, G. Hockncy, and A. Soni, Phys. Rev. D38, 3540 

(1988) . 

116. R. Sommcr, Phys. Rept. 275, 1 (1996) [hep-lat/9401037]. 

117. H. Wittig, Int. J. Mod. Phys. A12, 4477 (1997) [hcp-lat/9705034]. 

118. T. R. Klasscn, Nucl. Phys. B509, 391 (1998); Nucl. Phys. B Proc. Suppl. 
73, 918 (1999). 

119. J. Harada, A. S. Kronfcld, H. Matsufuru, N. Nakajima, and T. Onogi, 
Phys. Rev. D64, 074501 (2001) [hcp-lat/0103026]. 

120. S. Collins, et al, Phys. Rev. D64, 055002 (2001) [hcp-lat/0101019]. 



Lattice QCD 



65 



121. Z. Sroczynski, Nucl. Phys. B Proc. Suppl. 83, 971 (2000) [hep-lat/ 
9910004]; Int. J. Mod. Phys. A16S1C, 1231 (2001) [hcp-lat/0011059]. 

122. G. P. Lepage, Nucl. Phys. B Proc. Suppl. 26, 45 (1992). 

123. S. Hashimoto, Phys. Rev. D50, 4639 (1994) [hep-lat/9403028] . 

124. M. A. Shifman and M. B. Voloshin, Sov. J. Nucl. Phys. 45, 292 (1987); 
47, 511 (1988) [Yad. Fiz. 45, 463 (1987); 47, 801 (1988)]. 

125. N. Isgur and M. B. Wise, Phys. Lett. B232, 113 (1989); B237, 527 
(1990). 

126. P. A. M. Dirac, Proc. Roy. Soc. Lond. A133, 60 (1931). 

127. A. F. Falk, M. Neubert, and M. Luke, Nucl. Phys. B388, 363 (1992) 
[hep-ph/9204229]. 

128. Y. Kuramashi, Phys. Rev. D58, 034507 (1998) [hep-lat/9705036]. 

129. D. Becirevic et al, Phys. Rev. D60, 074501 (1999) [hep-lat/9811003]. 

130. K. C. Bowler et al. [UKQCD Collaboration], Phys. Lett. B486, 111 

(2000) [hep-lat/9911011]. 

131. K. C. Bowler et al. [UKQCD Collaboration], Nucl. Phys. B619, 507 

(2001) [hep-lat/0007020]. 

132. A. Abada et al, Nucl. Phys. B619, 565 (2001) [hcp-lat/0011065]. 

133. J. Rolf and S. Sint, Nucl. Phys. B Proc. Suppl. 106, 239 (2002) [hep- 
ph/0110139]. 

134. S. Hashimoto et al, Phys. Rev. D61, 014502 (2000) [hep-ph/9906376]. 

135. A. X. El-Khadra, A. S. Kronfcld, P. B. Mackenzie, S. M. Ryan, and 
J. N. Simone, Phys. Rev. D64, 014502 (2001) [hcp-ph/0101023]. 

136. J. Gasser and H. Lcutwyler, Ann. Phys. 158, 142 (1984). 

137. J. Gasser and H. Lcutwyler, Nucl. Phys. B250, 465 (1985). 

138. M. Gockeler et al, Phys. Rev. D62, 054504 (2000) [hep-lat/9908005]. 

139. A. Duncan, E. Eichten, and H. Thacker, Phys. Rev. Lett. 76, 3894 (1996) 
[hep-lat/9602005]. 

140. J. Garden, J. Hcitger, R. Sommer, and H. Wittig [Alpha Collaboration], 
Nucl. Phys. B571, 237 (2000) [hep-lat/9906013]. 

141. C. G. Boyd and B. Grinstein, Nucl. Phys. B442, 205 (1995) [hep- 
ph/9402340]. 

142. S. R. Sharpc and Y. Zhang, Phys. Rev. D53, 5125 (1996) [hep-lat/ 



9510037]. 

143. W. Bardecn, A. Duncan, E. Eichten, N. Isgur, and H. Thacker, Phys. 
Rev. D65, 014509 (2002) [hep-lat/0106008]. 

144. S. R. Sharpe, Phys. Rev. D41, 3233 (1990); 46, 3146 (1992) [hep- 
lat/9205020]. 

145. A. Morel, J. Phys. (France) 48, 1111 (1987). 

146. P. H. Damgaard, J. C. Osborn, D. Toublan, and J. J. Verbaarschot, Nucl. 
Phys. B547, 305 (1999) [hep-th/9811212]. 

147. S. R. Sharpe and N. Shoresh, Phys. Rev. D64, 114510 (2001) [hep- 



GG 



Handbook of QCD / Volume 4 



lat/0108003]. 

148. C. W. Bernard and M. F. L. Goltcrman, Phys. Rev. D46, 853 (1992) 
[hep-lat/9204007]. 

149. M. J. Booth, Phys. Rev. D51, 2338 (1995) [hcp-ph/9411433], hep- 
ph/9412228. 

150. S. R. Sharpe and N. Shoresh, Phys. Rev. D62, 094503 (2000) [hep- 
lat/0006017]. 

151. N. Kawamoto and J. Smit, Nucl. Phys. B192, 100 (1981). 

152. H. S. Sharatchandra, H. J. Thun, and P. Weisz, Nucl. Phys. B192, 205 
(1981). 

153. G. W. Kilcup and S. R. Sharpe, Nucl. Phys. B283, 493 (1987). 

154. H. Joos and M. Schafcr, Z. Phys. C34, 465 (1987). 

155. H. B. Nielsen and M. Ninomiya, Nucl. Phys. B185, 20 (1981); (E) 195, 
541 (1982). 

156. D. Fricdan, Commun. Math. Phys. 85 (1982) 481. 

157. P. H. Ginsparg and K. G. Wilson, Phys. Rev. D25, 2649 (1982). 

158. S. R. Sharpe and R. J. Singleton, Phys. Rev. D58, 074501 (1998) [hep- 
lat/9804028]. 

159. G. Rupak and N. Shoresh, hcp-lat/0201019. 

160. W. J. Lee and S. R. Sharpe, Phys. Rev. D60, 114503 (1999) [hep- 
lat/9905023]. 

161. C. W. Bernard and M. F. L. Goltcrman, Phys. Rev. D49, 486 (1994) 
[hep-lat/9306005]. 

162. K.-H. Mutter and K. Schilling, Nucl. Phys. B200, 362 (1982); Phys. Lett. 
B117, 75 (1982). 

163. M. Guagnclli, J. Hcitgcr, R. Sommcr, and H. Wittig [Alpha Collabora- 
tion], Nucl. Phys. B560, 465 (1999) [hep-lat/9903040] . 

164. M. Creutz, Phys. Rev. D15, 1128 (1977). 

165. M. Liischer, Commun. Math. Phys. 54, 283 (1977). 

166. K. Osterwalder and E. Seiler, Ann. Phys. 110, 440 (1978). 

167. H. Satz, Rept. Prog. Phys. 63, 1511 (2000) [hep-ph/0007069]. 

168. E. S. Fraga, R. D. Pisarski, and J. Schaffner-Bielich, Phys. Rev. D63, 
121702 (2001) [hep-ph/0101143]. 

169. S. A. Gottlieb et a/., Phys. Rev. Lett. 55, 1958 (1985). 

170. N. H. Christ and A. E. Terrano, Phys. Rev. Lett. 56, 111 (1986). 

171. R. V. Gavai et al. [MT C Collaboration], Phys. Lett. B232, 491 (1989); 
Phys. Lett. B241, 567 (1990). 

172. F. R. Brown et al, Phys. Lett. B251, 181 (1990); Phys. Rev. Lett. 65, 
2491 (1990). 

173. M. Fukugita, H. Mino, M. Okawa, and A. Ukawa, Phys. Rev. Lett. 65, 
816 (1990); Phys. Rev. D42, 2936 (1990). 

174. G. Boyd et al, Phys. Rev. Lett. 75, 4169 (1995) [hep-lat/9506025]; Nucl. 



Lattice QCD 



67 



Phys. B469, 419 (1996) [hep-lat/9602007]. 

175. H. Meyer-Ortmanns, Rev. Mod. Phys. 68, 473 (1996) [hcp-lat/9608098]; 
E. Lacrmann, Phys. Part. Nucl. 30, 304 (1999) [Fiz. Elem. Chast. Atom. 
Yadra 30, 720 (1999)]; S. Ejiri, Nucl. Phys. B Proc. Suppl. 94, 19 (2001) 
[hep-lat/0011006]; S. Hands, ibid. 106, 142 (2002) [hep-lat/0109034]. 

176. P. van Baal, in At the Frontier of Particle Physics: Handbook of QCD, 
Vol. 2, edited by M. Shifman (World Scientific, Singapore, 2001) [hep- 
ph/0008206]. 

177. M. Liischer, Commun. Math. Phys. 104, 177 (1986). 

178. M. Liischer, Commun. Math. Phys. 105, 153 (1986). 

179. M. Liischer, Nucl. Phys. B354, 531 (1991). 

180. M. Liischer, Nucl. Phys. B364, 237 (1991). 

181. J. Gasscr and H. Lcutwylcr, Nucl. Phys. B307, 763 (1988). 

182. P. Hasenfratz and H. Lcutwylcr, Nucl. Phys. B343, 241 (1990). 

183. F. C. Hansen, Nucl. Phys. B345, 685 (1990); F. C. Hansen and 
H. Lcutwyler, ibid. 350, 201 (1991). 

184. U. J. Wiese, Nucl. Phys. B Proc. Suppl. 9, 609 (1989). 

185. T. A. DcGrand, Phys. Rev. D43, 2296 (1991). 

186. L. Lellouch and M. Liischer, Commun. Math. Phys. 219, 31 (2001) [hep- 
lat/0003023]; See also Sec. 7 of L. Lellouch, Nucl. Phys. B Proc. Suppl. 
94, 142 (2001) [hep-lat/0011088], for a simple summary. 

187. L. Maiani and M. Testa, Phys. Lett. B245, 585 (1990). 

188. C. W. Bernard and M. F. L. Golterman, Phys. Rev. D53, 476 (1996). 

189. C. J. Lin, G. Martinclli, C. T. Sachrajda, and M. Testa, Nucl. Phys. 
B619, 467 (2001) [hep-lat/0104006]. 

190. A. Ali Khan et al. [CP-PACS Collaboration], Phys. Rev. D65, 054505 
(2002) [hep-lat/0105015]. 

191. J. Sexton, A. Vaccarino, and D. Weingarten, Phys. Rev. Lett. 75, 4563 
(1995) [hep-lat/9510022]. 

192. G. S. Bali, in Quark Confinement and the Hadron Spectrum IV, edited 
by W. Lucha and K. M. Maung (World Scientific, Singapore, 2001) [hep- 
ph/0010032]. 

193. C. J. Morningstar, K. J. Juge, and J. Kuti, Nucl. Phys. B Proc. Suppl. 
73, 590 (1999) [hep-lat/9809098]. 

194. J. Frohlich and P. A. Marchctti, Nucl. Phys. B Proc. Suppl. 106, 47 
(2002) [hcp-lat/01 10193]. 

195. D. E. Groom et al. [Particle Data Group], Eur. Phys. J. C15, 1 (2000). 

196. S. Chandrasekharan, Phys. Rev. D60, 074503 (1999) [hep-lat/9805015]. 



